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A.  Statement  of  the  Task  and  General  Approach 

This  report  documents  the  results  of  the  work  performed  under  the 
Pershing  II  Task  Order  DAAH01-81-D-AQ03,  Delivery  Order  0038  (GT/EES 
Project  No.  A-3165) .  The  task  order  required  that  the  following  subtasks 
be  performed: 

e  The  range  safety  problems  which  would  arise  if  the  motor  nozzle 
were  unintentionally  deflected  during  boost  phase  were  to  be  evaluated. 

e  The  aerodynamic  simulation  model  used  in  a  Pershing  II  computer- 
simulation  program  was  to  be  validated  against  actual  flight  test  data. 

•  Computer  simulation  results  were  to  be  compared  with  actual  flight 
data  for  a  set  of  Tactical  Ballistic  Missile  (TBM)  trajectory  profiles. 

The  results  of  these  comparisons  were  “to  be  presented  in  a  form  use¬ 
ful  for  evaluation". 

B.  Background 

The  task  documented  with  this  report  was  an  outgrowth  of  an  earlier 
task  (Pershing  II  Debris  Studies,  DAAH01-81-D-A003,  Delivery  Order  0017, 
GT/EES  Project  No.  A-2946).  The  earlier  task  was  aimed  at  evaluating 
launch  and  target  sites  for  Pershing  II  firings,  determining  the  effect  of 
varying  aerodynamic  conditions  on  the  missile's  flight  behavior,  and  con¬ 
ducting  range  safety  studies.  A  follow-on  task  order  request  is  in  pre¬ 
paration  and  will  expand  upon  the  work  performed  in  support  of  these  two 
projects. 

C.  Contents  of  This  Report 

A  detailed  discussion  of  the  subtask  objectives  and  accomplishments 
is  given  in  Section  II.  Section  III  discusses  the  U70  and  TRW  simulation 
programs  used  to  carry  out  Pershing  II  missile  simulations  [1].  It  in¬ 
cludes  a  description  of  the  U70  program  and  the  computer  on  which  it  runs, 
modifications  that  had  to  be  made  to  it  to  support  this  task,  and  the 
operational  procedures  for  running  it.  It  also  discusses  the  problems 
involved  in  converting  the  TRW  program  to  run  on  the  U.  S.  Army  Missile 
Command's  Perkln-Elmer  3220  computer.  Section  IV  discusses  the  results 
that  were  obtained  in  the  execution  of  this  task.  Finally  Section  V  pre¬ 
sents  recommendations  for  future  work. 


II.  DETAILED  OBJECTIVES 

A.  Subtask  Objectives  and  Accomplishments 

During  the  performance  of  this  task,  two  of  the  three  required  com¬ 
puter  simulation  studies  were  performed,  using  the  U70  trajectory  simula¬ 
tion  program  [1]  and  the  TKAJ  targeting  program.  The  third  study,  the 
comparison  of  Pershing  II  computer  simulation  studies  with  actual  flight 
data,  could  not  be  carried  out  because  flight  test  data  were  not  available 
during  the  period  of  performance  of  the  task.  The  effort  that  would  nor¬ 
mally  have  been  expended  upon  this  comparison  of  flight  test  data  and  simu¬ 
lation  results  was  directed  instead  toward  the  Implementation  of  a  Maximum 
Likelihood  Method  computer  subroutine  for  data  analysis. 

Similar  studies  were  to  be  undertaken  using  a  new  Pershing  II  simula¬ 
tion  program  developed  by  TRW,  Inc.,  but  the  full  computer  program  could  not 
be  obtained  In  time  for  use  on  this  task.  Its  treatment  in  this  report 
is  necessarily  restricted  to  a  discussion  of  the  conversion  and  docu¬ 
mentation  of  the  program  Itself  rather  than  the  results  obtained  using  it. 

A  more  detailed  description  of  the  four  subtask  objectives  is  given  below. 

1.  Flight  Safety  (Motor  Nozzle  Deflection)  Study 

The  motor  nozzle  deflection  study  was  a  flight  safety  study  con¬ 
ducted  to  determine  the  breakup  point  of  the  Pershing  II  missile  in  the 
event  of  an  engine  malfunction.  A  nozzle  deflection  would  cause  the 
missile  to  deviate  from  its  normal  flight  path.  The  objective  of  the  study 
was  to  determine  how  quickly  the  missile  would  leave  its  flight  corridor  in 
the  event  of  a  spurious  nozzle  deflection  and  whether  Range  Safety  would 
detect  this  aberration  in  time  to  permit  destruction  of  the  missile  before 
it  could  endanger  life  or  property. 

2.  Aerodynamic  Model  Validation  Study 

A  validation  of  the  aerodynamic  mathematical  model  used  in  the 
U70  Pershing  II  digital  simulator  was  required  for  utilization  of  flight 
data.  To  compare  the  simulation  results  to  the  flight  data,  a  statistical 
technique  called  the  Maximum  Likelihood  Method  was  to  be  used.  This  tech¬ 
nique  de-serves  special  mention  because  of  its  key  role  in  the  performance 
of  this  task.  This  technique  is  used  to  fit  theoretically  predicted  values 
of  relevant  variables  to  experimentally  derived  values  when  the  experimen¬ 
tal  values  contain  random  measurement  errors.  Generally,  the  Maximum 
Likelihood  Method  provides  better  estimates  of  the  underlying  true  values 
of  variables  than  a  least-squares  curve  fit  because,  in  making  its  estima¬ 
tes,  the  Maximum  Likelihood  Method  uses  a  prior  knowledge  of  the  standard 
deviations  of  instrument  measurement  errors  as  additional  information  to 
help  refine  its  estimates. 
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3.  Tactical  Ballistic  Missile  Trajectory  Study 

The  Tactical  Ballistic  Missile  (TBM)  study  took  the  form  of  a 
matrix  of  U70  runs  with  various  flight  paths  along  one  dimension  and  dif¬ 
ferent  flight  characteristics  along  the  other  dimension. 

One  launch  site  and  three  target  sites  were  used  to  generate  the 
three  flight  paths.  Then  for  each  flight  path,  a  set  of  three  runs  was 
made  for  three  assumed  conditions  -  first,  a  "nomln&x  motor"  or  unperturbed- 
baseline  flight  plan;  second,  a  ballistic  re-entry  vehicle  flight  plan  in 
which  no  aerodynamic  corrections  are  applied  to  the  missile  during  its  ter¬ 
minal  re-entry  phase;  and  third,  a  flight  plan  assuming  offsets  in  target 
latitude  and  longitude  from  its  normal  target  location.  Thus,  a  total  of 
nine  runs  were  made  for  this  study. 

B.  Conversion  and  Installation  of  The  TRW  Simulator 

The  TRW  program  is  a  six-degree-of-freedom  simulator  developed  for 
the  purpose  of  validating  Pershing  II  onboard  software.  This  simulator 
contains  aerodynamic  modeling  features  and  provides  simulation  capabilities 
not  found  in  the  U70  simulation  program  -  capabilities  which  the  U.  S.  Army 
Missile  Command's  Systems  Simulation  and  Development  Directorate,  Systems 
Evaluation  Branch,  seeks  to  add  to  its  repertoire  of  simulation  programs. 

The  job  of  understanding  the  program  and  converting  it  to  the  Missile 
Command's  computer  was  assigned  to  Georgia  Tech  as  a  part  of  this  task 
(A-3I65).  TRW  delivered  a  copy  of  its  program  to  the  Systems  Evaluation 
Branch.  However,  when  program  compilation  and  link  loading  was  attempted, 
certain  sections  of  the  program,  such  as  the  MAIM  control  program  and  the 
INTERP  routine  were  found  to  be  missing.  These  could  not  be  obtained 
during  the  period  of  performance  of  the  task,  and  therefore,  could  not  be 
converted.  Insofar  as  possible,  the  remainder  of  the  program  was  partially 
converted  from  TRW's  Digital  Equipment  Corporation  VAX  11/780  computer  to 
the  System  Simulation  and  Development  Directorate's  Perkin-Elmer  (PE)  3220 
computer  after  making  minor  changes  to  accommodate  for  the  differences  bet¬ 
ween  the  FORTRAN  capabilities  of  the  two  computers. 

Although  the  installation  of  the  TRW  program  on  the  PE  3220  computer 
and  the  incorporation  of  its  aerodynamic  model  within  the  1)70  program  is 
not  explicitly  required  in  this  task's  scope  of  work,  it  is  implicitly 
required  for  the  proper  performance  of  the  task.  Since  it  has  consumed 
much  of  the  task's  effort,  it  will  be  discussed  in  this  report  as  though  it 
were  a  fourth  task  objective. 
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III.  THE  U70  AND  TRW  SIMULATION  PROGRAMS 
A.  The  U70  Program 

1.  Description  of  the  U70  Program 

The  U70  trajectory  program  employed  in  the  performance  of  this 
task,  uses  the  given  prelaunch,  target,  and  flight  conditions  to  simulate 
the  flight  of  a  single-stage  or  two-stage  missile  and  its  associated 
maneuvering  re-entry  vehicle.  This  program  utilizes  the  equations  and 
requirements  found  in  the  Pershing  Launch  Computer  and  the  Pershing 
Airborne  Computer  to  simulate  the  missile  flight  from  launch  through  the 
boost,  midcourse,  and  terminal  portions  of  flight  to  impact  [2]. 

The  U70  trajectory  program  runs  on  a  Perkin-Elmer  3220  computer 
system  equipped  with  800  kilobytes  of  random  access  memory,  a  nine-track 
tape  drive,  three  67-megabyte  disk  drives,  and  five  demand  terminals 
(Including  the  system  console).  The  simulation  is  coded  in  FORTRAN  VII-D, 
although  it  is  probably  compatible  with  other  versions  of  FORTRAN. 

The  program  is  able  to  restart  the  missile  simulation  at  any 
time  during  the  terminal  re-entry  portion  of  flight  upon  entry  of  data  such 
as  time,  position,  velocity,  and  orientation  of  the  missile  at  the  desired 
restart  time,  plus  the  restart  data  acquired  from  a  previous  simulation. 

2.  Modifications  to  the  U70  Program 

a.  Modifications  for  the  Motor  Nozzle  Deflection  Study 

For  this  study,  modifications  were  made  to  the  U70  program's 
BAPLT  (Boost  Autopilot)  subroutine  to  enable  desired  deflection  angles  to 
be  input  at  some  time  during  boost.  These  modifications  are  as  follows: 

MODIFICATIONS  TO  BAPLT 

470  CONTINUE 

IPND  -  ICON(50) 

IYND  -  IC0N(51) 

IF( IPND  .EQ.  0)  GO  TO  453 
IF  (TIMC  .LT.  CIN(382))  GO  TO  453 
DNV  (1,  1)  -  CIN  (380)  *  DTOR 

453  CONTINUE 

IF( IYND  .EQ.  0)  GO  TO  454 
IF(TIMC  .LT.  CIN(382))  GO  TO  454 
DNV  (1,  2)  -  CIN(381)  *  DTOR 

454  CONTINUE 
RETURN 
END 


4 


b.  Aerodynamic  Model  Flight  Validation 


As  mentioned  previously,  the  Maximum  Likelihood  Method  (MLM) 
of  data  correlation  was  intended  to  help  fit  the  simulator-generated  aero¬ 
dynamic  flight  data  to  the  actual  test-flight  telemetry  data,  i.e.,  to 
"average  out"  random  telemetry  Instrumentation  errors  in  order  to  provide 
more  meaningful  comparisons  between  the  simulated  and  experimental  results. 
The  machinery  for  accomplishing  this  was  incorporated  into  a  new  U70 
subroutine  called  DATACORR.  The  addition  of  DATACORR  to  the  070  program 
required  the  creation  of  a  new  link  procedure,  TU70LINK.CSS,  which  builds  a 
new  task  called  TU70. 

The  TU70  task  requires  as  inputs  the  instantaneous  inertial- 
frame  X(t),  Y(t),  and  Z(t)  position  coordinates  and  $(t),  6(t),  and  <Kt) 
(roll,  pitch,  and  yaw)  angular  coordinates  of  the  missile  and  their  time 
derivatives  (where  t  represents  time).  For  program  checkout  purposes, 
constant  offsets  have  been  programmed  into  the  inertial  position  and  rota¬ 
tion  inputs.  Since  the  MLM  uses  an  iterative  technique,  an  upper  limit 
must  be  entered  on  the  number  of  iterations  that  the  computer  is  to  attempt 
before  giving  up  (in  case  the  iteration  sequence  does  not  converge). 

Finally,  a  6  x  6  matrix  of  standard  deviations  and  cross-correlations  of 
measurement  noise  must  be  entered  for  the  three  position  and  the  three 
rotation  coordinates  at  t  ■  0,  e.g.,  at  launch  (Table  1). 

This  matrix  contains  in  its  [1,  1]  location,  the  variance  (o^) 
of  the  instrumentation  measurement  errors  of  the  missile's  x-coordinate  at 
t  -  0.  In  its  [1,  2]  location,  it  contains  the  mean  square  error  in  X  due 
to  an  error  in  Y  and  so  forth,  for  all  36  elements.  These  TU70 

input  requirements  are  summarized  in  Table  2.  The  coding  for  the  DATACORR 
subroutine  is  reproduced  in  Appendix  A. 

c.  Data  Format  Modifications  for  the  Tactical  Ballistic 
Missile  Simulations 

Modifications  wepe  made  for  the  Tactical  Ballistic  Missile 
trajectory  simulations  to  provide  the  same  output  format  for  all  stages  of 
flight.  Different  variables  were  output  for  the  boost  and  for  the  re-entry 
stages;  these  variables  were  combined  to  produce  one  block  of  output  data 
for  both  stages.  This  required  changes  in  the  BOUT  (Boost  Output)  and 
RVOUT  (Re-entry  Vehicle  Output)  subroutines.  These  changes  are  documented 
in  Appendix  B. 
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TABLE  1.  COVARIANCE  MATRIX  OF  MEASUREMENT  NOISE,  N 
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TABLE  2.  TU70  INPUTS  FOR  THE  DATA  CORRELATION  FEATURE 


INDICATORS: 

CIN(75)  -  Contains  an  on-off  "switch”  for  the  data  correlation  feature 
(0  •  Off,  1  ■  On) 

CIN(70)  -  Contains  an  on-off  "switch”  for  the  computation  of  the  transpose 
of  the  sensitivity  matrix,  S  (0  ■  Off,  1  -  On) 


CONSTANTS: 
CIN(499)  - 
CIN(498)  - 
CIN(497)  - 

CIN(496)  - 
CIN(495)  - 
CIN(494)  - 
CIN(493)  - 

CIN(492)  - 
CIN(49I)  - 
CIN(490)  - 
CIN(489)  ■ 

CIN(488)  - 
CIN<4 87)  - 


Contains  the  number  of  experimental  data  points 


°x2 


°y2 


I  Contains  the  initial  inputs  of  the  co- 
)  variance  matrix  of  the  measurement 
{  noise  (assumes  a  diagonal  matrix) 


Ee 

E* 

Eij; 

Ex 


Contains  the  error-vector  differences 
(  between  the  measured  and  computed  values 
of  x,  y,  z,  9,  $,  and  <l> 


3.  Instructions  for  Using  the  U70  Program 

A  description  of  the  procedures  for  compiling,  running,  linking, 
restarting  the  program,  and  making  output  tapes  is  given  in  Appendix  C, 
together  with  a  list  of  the  .CSS  files  used  to  accomplish  these  tasks,  and 
a  list  of  the  U70  source  files. 

The  TU70  program  contains  an  array  of  input  constants  and  indica¬ 
tors  called  CIN.  Table  2  contains  definitions  for  the  constants  and  indi¬ 
cators  which  serve  as  inputs  for  the  DATACORR  subroutines  for  the  TU70 
task. 


B.  The  TRW  Program 

1.  Compiling  the  TRW  Simulation  Program 

The  TRW  program  was  originally  set  up  for  interactive  com¬ 
pilation,  program  linkage,  and  execution.  However,  in  converting  the 
program  to  the  Perkln-Elmer  3220  computer,  it  became  necessary  to  switch  to 
batch  mode  compilation  and  linkage.  The  reason  for  turning  to  batch  mode 
operation  was  that  not  every  file  or  function  required  by  the  program  is 
present  in  the  program  itself.  Some  of  the  necessary  information  is  stored 
in  independent  files,  and  operating  in  batch  mode  permits  these  independent 
files  to  be  found  during  compilation  and  to  be  linked  to  the  MAIN  program 
during  the  program  linkage  phase.  To  accomplish  this  result,  the  $ BATCH 
command  had  to  be  placed  at  the  beginning  of  the  program,  the  $BEND  command 
at  the  end  of  the  program,  and  the  $PR0G  declarations  put  at  the  beginning 
of  each  subroutine. 

Next,  a  "compile  file"  was  produced  by  modifying  a  copy  of  the 
systems  file  F7CAE.CSS  (which  contains  the  FORTRAN  VII  compiler  and  the 
linking  loader)  and  then  using  it  to  compile  the  TRW  program.  Compilation 
time  was  saved  by  compiling  one  subroutine  at  a  time,  producing  an  object 
code  image  and  saving  it  in  a  separate  file,  and  then  linking  it  together 
with  other  object  code  subroutines  during  the  link  phase.  This  meant  that 
only  those  subroutines  which  had  been  updated  had  to  be  recompiled. 

F7CAE.CSS  is  a  procedure  which  compiles  and  links  any  program. 

The  compiler  in  F7CAE.CSS  uses  system  file  F7D.TSK  to  produce  an  object 
file  and  F7D.ERR  to  record  errors  found  during  the  compile  phase.  The 
F7CAE.CSS  link  loader  processes  object  files  generated  by  the  compiler  and 
creates  a  task  file  for  execution.  The  link  loader  also  creates  a  load  map 
and  provides  a  record  of  any  undefined  external  references  or  symbols. 
Appendix  D  lists  the  file  names  used  by  the  TRW  program  and  Appendix  E 
lists  the  program's  subroutine  names  and  calling  sequences. 
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Of  Che  many  error  messages  which  surfaced  during  the  first  com¬ 
pilation,  a  substantial  fraction  was  associated  with  VAX-peculiar  FORTRAN 
INCLUDE  statements.  INCLUDE  is  a  VAX  FORTRAN  compile-time  comtand  which 
causes  the  compiler  to  copy  a  disk  file  into  a  program  -  typically,  a  file 
containing  COMMON  and  EQUIVALENCE  statements.  This  permits  a  simple 
INCLUDE  statement  to  be  substituted  In  subroutines  for  the  extensive  COMON 
and  EQUIVALENCE  statements  that  would  ordinarily  be  found  there  instead. 
This  practice  not  only  reduces  the  amount  of  effort  required  to  write  and 
update  the  program  but  also  reduces  the  chances  of  making  a  mistake  in 
coding,  since  the  COMMON  and  EQUIVALENCE  statements  only  have  to  be  updated 
at  one  place  in  the  program  rather  than  in  every  subroutine  in  which  they 
are  referenced. 

Another  VAX-peculiar  FORTRAN  VII  enhancement  is  the  DO  WHILE  com¬ 
mand.  Since  the  DO  WHILE  statement  is  not  recognized  by  Perkin-Elmer's 
FORTRAN  VII  compiler,  it  was  necessary  to  replace  DO  WHILE'S  in  the 
original  program  with  ordinary  DO  statements  in  the  Perkln-Elmer  version  of 
the  program. 

The  above  modifications  constitute  the  principal  programing 
changes  that  have  so  far  had  to  be  made  to  the  TRW  subroutines  in  order  to 
convert  them  from  the  VAX  11/780  computer  to  the  Perkin-Elmer  3220  com¬ 
puter. 

2.  Subroutines  Missing  From  the  TRW  Program 

As  mentioned  in  Section  II. E.,  the  TRW  simulation  program  was 
delivered  with  some  major  portions  of  the  program  missing.  The  most  cri¬ 
tical  missing  subroutine  was  INTERPC Interpolation)  which,  by  interpolating 
large  data  tables,  would  have  provided  the  thrust  characteristics,  mass 
properties,  and  aerodynamic  data  needed  to  simulate  the  flight.  The  main 
control  program,  MAIN,  was  also  missing,  although  a  program  SDCTRL 
( Six-Degree-of-Freedom  Control)  was  found  which  seems  to  perform  similar 
functions.  Essential  portions  of  the  program  that  were  required  to  simu¬ 
late  steering,  guidance,  and  navigation  in  the  Pershing  Airborne  Computer 
were  also  missing.  Finally,  an  error  message  routine  ERRMSG  was  missing. 

a.  INTERPC Interpolation) 

With  the  aid  of  a  former  developer  of  the  TRW  INTERP  routine, 
a  prior  version  of  INTERP  was  located  in  an  earlier  edition  of  the  TRW 
simulator  and  efforts  were  made  to  understand  it  and  use  it  in  the  current 
version  of  the  TRW  program.  At  the  heart  of  the  INTERP  routine  is  the 
above  mentioned  Involved  set  of  aerodynamic  tables,  organized  in  a  way  that 
minimizes  storage  requirements  without  unduly  Increasing  run  times.  A 
discussion  of  what  has  been  learned  about  this  multi-dimensional  data  table 
is  presented  in  Appendix  F. 


b.  MAIM  and  SDCTRL  ( Six-Degree-of-Freedom-Control)  Programs 

After  finding  that  the  MAIN  routine  was  missing  from  the  TRW 
simulator,  efforts  were  begun  to  develop  such  a  function-controlling 
routine.  An  earlier  version  of  MAIN  was  found  in  a  program  listing  and  the 
listing  was  used  to  help  understand  what  was  needed  to  recreate  a  current 
MAIN  routine.  A  flowchart  of  this  early  version  MAIN  routine  is  contained 
in  Appendix  G.  Appendix  H  contains  a  flowchart  of  the  SDCTRL  routine. 

c.  The  Output  Processor 

The  TRW  simulation  program  has  an  output  processor  which  takes 
advantage  of  large  arrays  to  store  values  of  variables  for  output.  These 
arrays  are  saved  in  three  groups  named  BIG0LR,  BIG01I,  and  BIG01D,  which 
contain  Real,  Integer,  and  Double  Precision  variables,  respectively.  These 
output-variable  values  are  transferred  to  the  BIG01  arrays  for  use  by  the 
output  processor  through  EQUIVALENCE  statments  located  in  each 
subroutine  [1].  Each  variable  stored  in  BIG01R,  BIG01I,  and  BIG01D  is 
stored  on  a  disk  file  BIG01.DAT.  The  names  of  these  variables  are  listed 
in  Appendix  I. 

3.  Tape-to-Disk  Conversion  of  the  TRW  Program 

TRW  delivered  its  simulation  program  to  the  System  Evaluation 
Branch  on  a  nine-track  1600  bit-per-inch  tape.  This  was  loaded  onto  a 
Perkln-Elmer  3220  disk  pack  using  the  following  procedure: 

Mount  the  tape  on  a  1600  bpi  tape  drive. 

>  COPY 

>  AL  FILE1 .FTN,IN,80 

>  OUT  FILE1.FTN 

>  COPY  *,* 
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IV.  RESULTS 

A.  Flight  Safety  (Motor  Nozzle  Deflection)  Results 

As  mentioned  earlier  in  Section  II, B,  this  Pershing  It  missile  simu¬ 
lation  was  carried  out  using  the  U70  simulator  to  determine  what  would  hap¬ 
pen  if,  through  some  hardware  failure,  the  nozzle  were  accidently  deflected 
during  a  live  missile  firing.  In  carrying  out  these  simulated  flight 
tests,  the  missile  was  allowed  to  "fly"  unperturbed  (in  the  computer)  for  a 
short  time.  Then  the  motor  nozzle  was  deflected  by  a  pre-determined  angle 
and  the  simulation  continued  until  either  the  missile's  total  angle  of 
attack  exceeded  15°  or  the  normal  (perpendicular  to  the  body  of  the 
missile)  acceleration  exceeded  5  g's.  Either  one  of  these  occurences  was 
assumed  to  generate  sufficiently  unstable  conditions  that  breakup  of  the 
missile  would  occur,  so  the  simulation  was  terminated  at  that  point. 

The  spurious  nozzle  deflections  were  assumed  to  occur  at  30  seconds  and 
at  49  seconds  into  the  flight.  The  30-second  flight  time  was  chosen 
because,  at  30  seconds,  the  aerodynamic  forces  acting  on  the  missile  would  be 
at  or  near  their  maximum,  i.e,  a  worst-case  condition.  The  49-second 
flight  time  was  chosen  because  it  is  almost  at  first-stage  burnout,  and  is 
a  time  when  missile  failure  is  likely.  Three  nozzle  deflection  angles  were 
tried,  0.5°,  2°,  and  7.6°,  first  in  pitch  and  then  in  yaw,  leading  to  a 
total  of  twelve  cases  (two  deflection  times  and  three  nozzle  deflection 
angles,  first  in  pitch  and  then  in  yaw).  All  the  cases  were  simulated 
flights  from  McGregor,  New  Mexico,  to  McDonald’s  Well,  New  Mexico,  with  the 
targeting  conditions  given  in  Table  3  below.  Before  presenting  these 
results  in  detail,  some  data  comparison  problems  need  to  be  discussed. 

TABLE  3.  LAUNCH  SITE  AND  TARGET 


McGregor,  NM  LAUNCH  SITE 

Latitude  -  32.09575° 
Longitude  ■  -106.2035° 
Altitude  ■  1251.0000  m 


McDonald's  Well,  NM  TARGET  SITE 

Latitude  -  33.113580° 
Longitude  ■  -106.35897° 
Altitude  ■  1228.00000  m 
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In  order  to  provide  a  standard  by  which  to  compare  the  deviated- 
nozzle  runs,  a  nominal  or  baseline  simulation  was  run  in  which  it  was 
assumed  that  the  missile  followed  a  nominal,  undeflected  flight  plan. 

Since  the  deviated  nozzle  results  were  printed  at  0.1  second  Intervals,  it 
would  have  been  advantageous  to  have  also  printed  the  baseline  results  at 
0.1  second  intervals  for  purposes  of  comparison.  Unfortunately,  the  large 
volume  of  printout  generated  by  the  U70  program  made  it  impractical  to 
print  results  at  0.1  second  intervals,  and  the  results  were  available  only 
at  1.0  second  Intervals.  Consequently,  the  baseline  results  had  to  be 
interpolated  at  0.1  second  intervals  in  order  to  compare  them  to  the 
devlated-nozzle  results,  i.e.,  a  direct  comparison  was  not  possible. 

A  second  data  comparison  problem  resides  in  the  fact  that,  for 
unknown  reasons,  the  first  data  points  in  the  deviated-nozzle  simulation 
printouts  (the  results  for  30.0  seconds  and  the  results  for  49.0  seconds) 
are  invalid.  These  30-  and  49-second  numbers  are  not  random  but  seem  to 
come  from  some  earlier  time  in  the  flight.  Subsequent  results  appear  to  be 
correct  when  compared  to  the  unperturbed  flight  results.  For  example,  the 
interpolated  unperturbed  flight  results  for  30.1  seconds  agree  to  four 
decimal  places  with  the  deviated-nozzle  results,  with  subsequent  data 
departing  from  the  baseline  data  in  a  regular  and  reasonable  way.  It  is 
the  authors'  conclusion  that  the  faulty  first  point  in  each  set  of  numbers 
was  the  result  of  a  print-routine  malfunction  rather  that  an  error  in  the 
deviated  nozzle  results.  However,  the  reader  needs  to  be  aware  of  this 
characteristic  in  the  data. 

Figure  1  depicts  the  missile's  flight  path  for  an  unperturbed  flight 
plan  and  for  the  three  different  yaw-axis  nozzle  deflection  angles  looking 
down  from  above  when  the  nozzle  deflections  occur  30.0  seconds  into  the 
flight.  Figure  1  is  a  plot  of  downrange  coordinates  versus  crossrange 
coordinates.  The  "straight  line”  in  Figure  1  represents  the  path  of  an 
unperturbed  missile.  (It  curves  about  1°  to  the  right  but  the  curvature  is 
too  small  to  be  visible  in  the  plot.)  For  all  three  deflections,  the 
missile  disintegrates  because  the  yaw  angle  of  attack  exceeds  15°  rather 
than  because  the  normal  acceleration  exceeds  5  g's.  Note  that  for  the 
smaller  the  nozzle  deflection  angle,  the  farther  the  missile  can  travel 
before  it  reaches  the  critical  angle  of  15s. 

Figure  2  shows  a  similar  plot  when  the  nozzle  deflection  occurs  49 
seconds  into  the  flight.  Note  that  the  flight  path  deflections  appear  to 
be  smaller  than  they  were  when  the  nozzle  deflection  occurred  at  30  seconds 
(Figure  2).  In  reality,  the  deflections  are  the  same  but  the  velocity  of 
the  missile  is  greater,  and  this  increases  the  horizontal  scale  of  the 
plot.  (The  results  for  a  nozzle  deflection  of  7.6s  are  missing  because  of 
a  faulty  computer  run.) 
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Figure  1.  Crossrange  vs.  Downrange  for  Yaw  Nozzle  Deflections  (6's)  Occurring 
at  30  Seconds  Into  the  Flight. 
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DOWN RANGE  (m) 

Figure  2.  Crossrange  vs.  Downrange  for  Yaw  Nozzle  Deflections  (6's)  Occurring  at 
49  Seconds  into  the  Flight. 


Figure  3  shows  a  plot  of  Che  missile  vertical  pitch  angle  versus  its 
downrange  coordinates  when  a  nozzle  deflection  in  pitch  occurs  at  30 
seconds,  and  Figure  4  depicts  similar  plots  for  the  49  second  pitch  nozzle 
deflection. 

Figures  S  through  8  show  how  the  normal  (lateral)  accelerations  of 
the  missile  vary  with  time.  An  examination  of  these  curves  confirms  that 
in  every  case,  the  simulated  flights  terminated  before  the  normal  accelera¬ 
tions  exceeded  the  destructive  limit  of  5  g's,  i.e.,  the  flights  terminated 
because  the  missile's  total  angle  of  attack  exceeded  15a.  There  is  an 
interesting  tendency  for  these  normal  accelerations  to  rise,  dip,  and  then 
rise  again.  No  simple  explanation  has  been  advanced  for  this  behavior, 
although  the  Pershing  II  guidance  system  is  sufficiently  complex  that  a 
complex  response  might  be  expected  in  the  event  of  a  major  nozzle  deflec¬ 
tion  malfunction.  (The  results  for  a  nozzle  deflection  of  7.6°  are  missing 
because  of  the  faulty  computer  run  mentioned  previously  in  connection  with 
Figure  2.) 

These  simulations  show  that  the  missile  would  leave  its  allowable 
flight  corridor  and  exceed  the  allowable  angle  of  attack  within  approxi¬ 
mately  one  second. 

B.  Aerodynamic  Model  Validation  Results 

As  mentioned  in  Section  II. B.,  this  study  could  not  be  completed 
because  Pershing  II  flight  test  data  were  not  available  in  time. 

In  support  of  this  task,  the  DATACORR  subroutine,  incorporating  the 
Maximum  Likelihood  Method  of  data  correlation,  was  written,  compiled,  and, 
insofar  as  possible,  was  tested.  However,  one  important  part  of  this 
subroutine,  the  derivation  of  sensitivity  coefficients,  was  not  completed. 

C.  Tactical  Ballistic  Missile  Trajectory  Results 

As  mentioned  in  Section  II. C. ,  nine  U70  runs  were  made  for  this  study 
for  nine  input  trajectory  profiles.  Three  of  these  profiles  entailed  off¬ 
sets  in  target  latitude  and  longitude. 

When  the  study  was  completed,  the  results  were  output  on  tape  and  deli¬ 
vered  to  the  sponsor. 
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DOWN RANGE  (m) 

Figure  3.  Pitch  Body  Angle  vs.  Downrange  for  Pitch  Nozzle  Deflections  (S' a)  Occurring 
at  30  Seconds  into  the  Flight. 
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DOWN RANGE  (m) 

Figure  4.  Fitch  Body  Angle  ve.  Dovnrange  for  Fitch  Nozzle  Deflections  (6's)  Occurring 
at  49  Seconds  into  the  Flight. 
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Figure  7.  Normal  Accelerations  vs.  Time  for  Yaw  Nozzle  Deflections  (6's)  Occurring 
at  30  Seconds  Into  the  Flight. 
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Figure  8.  Normal  Accelerations  vs.  Time  for  Yaw  Nozzle  Deflections  (6's)  Occurring 
at  49  Seconds  into  the  Flight. 


D.  Conversion  and  Installation  of  the  TRW  Simulation  Program 

Some  of  Che  deCails  of  the  TRW  program  conversion  have  been  discussed 
in  Section  III.B.,  and  in  Appendices  D  through  I.  A  summary  of  the  results 
is  as  follows: 

1.  All  subroutines  except  those  missing  from  the  incoming 
program  tape  were  converted  and  successfully  compiled. 

2.  The  organization  of  the  INTERP  tables  was  deciphered  and 
documented  (Appendix  F). 

3.  The  MAIN  and  SDCTRL  subroutines  were  flowcharted  (Appendices 

G  and  H). 

4.  A  linking  task  builder  TRWLINK.CSS  and  a  compiler  task  builder 
DAVEF7.CSS  were  created  to  accommodate  the  TRW  program. 
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V.  RECOMMENDATIONS  FOR  FURTHER  WORK 

A.  Flight  Safety  (Motor  Nozzle  Deflection)  Recommendations 

One  of  the  recommendations  for  future  work  is  to  carry  out  these 
evaluations,  re-running  the  simulation  with  print  statements  that  monitor 
the  moment-by-moment  variations  in  the  angles  of  attack,  the  down-range  and 
cross-range  coordinates,  and  the  inertial  flight  path  angles  before, 
during,  and  after  the  moment  when  the  nozzle  deviation  occurs. 

B.  Aerodynamic  Model  Validation  Recommendations 

Recommendations  for  future  work  consist  of  completing  the  DATACOR& 
subroutine  and  carrying  out  the  study  using  actual  flight  data. 

Normally,  in  employing  the  Maximum  Likelihood  Method,  all  the  com¬ 
puted  constitutive  forces  and  moments  that  make  up  the  total  X,  Y,  or  Z 
forces  or  <)> ,  0 ,  or  i>  moments  (constituents  of  lift,  drag,  etc.,)  are  com¬ 
pared  individually,  with  their  equivalent  measured  values.  However,  the 
aerodynamic  model  used  in  the  U70  simulator  is  not  yet  sufficiently 
detailed  to  permit  making  these  comparisons  because  of  the  difficulty  of 
attributing  corrections  to  individual  components  of  the  U70  predictions,  so 
the  comparisons  are  to  be  made  at  the  level  of  overall  body  forces  and 
torques.  For  this  reason,  the  DATACORR  subroutine  doesn't  presently  output 
Its  force  and  torque  corrections  to  the  main  U70  program.  However,  future 
plans  call  for  the  upgrading  of  the  main  program  to  provide  aerodynamic 
modeling  at  the  constitutive  level.  When  that  Is  done,  the  DATACORR 
subroutine  and  the  TU70  task  will  be  upgraded  to  output  these  corrections 
to  the  U70  program. 

C.  Tactical  Ballistic  Missile  Trajectory  Recommendations 

Recommendations  for  further  work  depend  upon  the  sponsor's  require¬ 
ments. 


D.  Conversion  and  Installation  of  the  TRW  Simulation  Program 

The  Maximum  Likelihood  Method  correlation  technique  being  used  to 
validate  the  U70  simulator  predictions  against  actual  flight  data  (Section 
II. C.)  should  be  used  to  evaluate  the  TRW  simulator. 

Recommendations  for  future  work  consist  of  completing  the  conversion 
and  Installation  of  the  TRW  program,  using  it  to  simulate  Pershing  II  test 
flights,  and  comparing  its  results  with  the  actual  test  flight  data. 
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APPENDIX  A 

CODE  LISTING  FOR  THE  DATACORR  (Data  Correlation)  SUBROUTINE 
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no  -nJ  r>  v*  uy  rw 


TABLE  A-l.  CODE  LISTING  FOR  THE  DATACORR  SUBROUTINE 

5fros  !*0mrc7P  - 

♦D3R0UTINB  BAT ACORR 
UNCLLDE  9/U70TIM.COM 
Z  ***  9EGIN  COMMON  PACK  *** 

IMPLICIT  D0U3LE  PRECISION  (A-H/C-Z) 

REAL  C3V/CUR3V/DCB/DCI/SCD/SCI/PSAVE 


COMMON 

C IN (400), SCO (HOC)/ SCI C11 00) 

/ DC9C7200) /DC 

I (&40) 

COMMON 

TITL(12)/TMAJ(3Q)/TAB(50)/TTA3(9C) 

/ CU  R3V  < 

3300) / DUM2  (1  0) 

COMMON 

TCCN(10Q)/ID6(32)/IDC(33)/ID 

1(32)/ 

ISC  (31  ) 

/ISI(30)/IND3V(1 

COMMON 

ND9(32)/MDI(3 

2)/NSI(30)/N3TR 

(30)/NTV(3/10)/IALISN 

COMMON 

/CONST/  A13 

/A23 

/  A33 

/  C  0  SB  L 

/COSPL 

/ CPSIA  / CTHA 

0 

CXNUP 

/CXNUY 

/DTCR 

/ 02AC2 

/  F1 1 

/  Fl 2  / FI  3 

0 

F  21 

/e22 

/  F  23 

/  F31 

/  F  32 

/ F3  3  /RL 

0 

R»RES 

/  RTOD 

/ S IN9L 

/ S INPL 

/SPSIA 

/STHA  / SXNUP 

0 \ 

SXNUY 

^  S1 1 

/Si  2 

/  Si  3 

/S21 

/S22  /S23 

0* 

S31 

/S32 

/  S33 

/TANPL 

/TOLH 

/T0L1D  /T0L2D 

0 

U11 

/  U2 1 

/  U22 

/  U2  3 

/U31 

/U33  /PI 

0 

COT 

/SOT 

/PIC2 

/TWOPI 

/THIRD 

/ S8RT2  / SQRT3 

0 

P 104 


COMMON 

AC  A 

/ACCEL 

/ACR 

/  ACR4 

/ACR5 

/ 3PINT 

/BYINT 

/  CGDE 

COMMON 

D3AMS 

/DG3T 

/ORHO 

/OVEL 

/ G AMMC 

/3AMWP 

/GAMRE 

/IG'JID 

COMMON 

I  -*0  0 

/ITHM 

/.I  YAW 

/PHIPE  /PHIRE 

/PHIYE 

/PHPEC 

/ PHPEP 

COMMON 

FHYEC 

/PHYE° 

/PSIC 

/PX 

/PY 

/®Z 

/  RCX 

/  RC  Y 

common 

RCZ 

/RD 

/ROD 

/  RM  AG 

/SGDE 

/SGE 

/  S  3G 

/TFF 

COMMON 

TFFS 

/  TG  E 

/TGN 

/  TGS 

/THETA 

C/TIMPL 

/URANX 

/URANY 

COMMON 

URANZ 

/UR  AX 

/URAY 

/  JRAZ 

/URAZX 

/ JR  AZ Y 

/URAZZ 

/  URX 

COMMON 

LRY 

/UPZ 

/VE 

/VEX 

/VEXS 

/  VE  Y 

/VEYS 

/VE  Z 

common 

VEZ3 

/VI® 

/  VLEX 

/  VL  E  Y 

/VLEZ 

/VR 

/VRE 

/  VREH 

COMMON 

COMMON 

V  °E  HX 

Z  INTS 

/VREHY 

/ V  RE  HZ 

/ ZDIF  F  /ZETAE 

/ZETR5 

/ Z INTC 

/ZINTP 

COMMON 

A  A1 

/  AA2 

/ALDDTC 

/ ALDDTS/ALOAD 

F/ALOADS 

/ALPDOE/ALPHDD 

COMMON 

ALPHP 

/ AL®HPD 

/ ALPHPS 

/ALPHS-  / ALPHT 

/ALPHY 

/ ALPHYO 

/ALPHYS 

COMMON 

ATC 

/ATP 

/AX 

/AY 

/AZ 

/All 

/A12 

/  A21 

COMMON 

A  22 

/A31 

/  A  32 

/  3B  2 

/  3  E  0  OT 

C/3E0DTS 

/BETADO 

/3ETAP 

COMMON 

ccTAR 

/3ETAY 

/ 3  ET  AW 

/3E7DDE/BTC 

/  5T  P 

/B11 

/  31  2 

common 

=  15 

/9?1 

/S22 

/  3  2  3 

/  i  31 

/  33  2 

/B3  3 

/CA 

COMMON 

C  A  A  ° 

/  C  A  A  Y 

/CAIN 

/CLP 

/  C  MQ 

/  C  M  XP 

/CN 

/CNAPO 

COMMON 

CNAYD 

/  CN  P 

/  C  NC 

/ CNVDO  /CNVDP 

/  CN  Y 

/  CNZ 

/ COSPSD 

COMMON 

C  OS  Trl  D 

/CPHI 

/CPSI 

/  CTHE 

/CY 

/CY  P 

/CYR 

/C11 

COMMON 

C  1  2 

/ Cl  3 

/DELTA1 

/0ELTA2/DELTA 

3/0ELTA4 

/  D  FP 

/OFR 

COMMON 

D  F  Y 

/OT 

/  D  VT  C 

/  OVTP 

/  0  VT  S 

/Dll 

/  D1  2 

/  Dl  3 

COMMON 

C  21 

/  02 : 

/D  23 

/  031 

/  D  32 

/  03  3 

/  DELC 

/  EDLC 

COMMON 

tODIC 

/-PMDC 

/EPHOS 

/EPHIC  /EPHIP 

/EPHIS 

/EPSDC 

/ EPSDS 

common 

E  °S  I C 

/EPSI° 

/cPSIS 

/ETA 

/ETHAC 

/ETHAP 

/ETHAS 

/ ETHD  C 

common 

ETHOS 

/FS 

/  FP 

/FROLL  / F  SL 

/FT 

/FX 

/  FXM 

COMMON 

F  Y 

/  F  V  A  V 

/  F  Y  A  w 

/  CYM 

/F  Z 

/FZAV 

/f:m 

/G 

COMMON 

G  A 

/GAODTC 

/GACOTS 

/GAMDO  /GAMOC 

E/3R 

/GTC 

/  GTP 

COMMON 

H 

/MOT 

/  H 1 1 

/Hi  2 

/HI? 

/  H21 

/  H  22 

/  H23 

COMMON 

M  51 

/  H?2 

/H*3 

zIAFRO  / I CN 

/ICO 

/ 1 0 1  NT 

/IDISC 
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TABLE  A-l. 

COMMON  1 0 T 
COMMON  IHI J 
COMMON  I®LOT 
COMMON  ITA 
COMMON  IUPD 
COMMON  LINC 
COMMON  Pt>$ 
COMMON  PS 
COMMON  PYPC 
COMMON  P 1 3 
COMMON  «  A 
COMMON  R  AIM 
COMMON  sigma 
COMMON  T 30 
COMMON  7  MET  A 
COMMON  TSTMA 

common  vtyp 

COMMON  VPRIY 
COMMON  V VIE 
COMMON  WD 
COMMON  wpes 
COMMON  XDG 
COMMON  XIYY 
COMMON  X LI  2 
COMMON  X  M 

common  xmavz 

COMMON  XMTZ 
COMMON  X  NV  3 
COMMON  X  3 
COMMON  YDS 
COMMON  Y  °D  S 
COMMON  ZDG 
COMMON  Z? 
COMMON  ZZ1 

COMMON  ACC 
COMMON  ADE32 
COMMON  AlYC 
COMMON  APX 
COMMON  CAeAS 
COMMON  CILYA 
COMMON  DVD 
COMMON  GADYC 
COMMON  GDPC 
COMMON  GLIM 
COMMON  IOUM 
COMMON  I  SEA 
COMMON  PHTTA 
COMMON  fi  <E 
COMMON  TLAM 
COMMON  V N A V 


CODE  LISTING  FOR  THE  DATACOSR  SUBROUTINE  (Continued) 


/ ICUSL  /I ENV 

/ IFAIL 

+ZFZO* 

/ I F INC 

/I FLAG 

TTGC — 

/IHT  /IMAJ 

/IMISS 

/INCH 

/INDT 

/INT 

/IPAGE 

/ IPNT  / 1 PP  NT 

/ IP  RC  3 

/  I  0  J  u 

/IRUN 

/  ISEN 

/ISNTH 

/  IT  3  /  ITC 

/ITHP 

/ IT0L8 

/ ITPOS 

/ITT 

/IUHT 

/IVAC  /IWD 

/  IW  IN 

/  J  STS 

/<FLAG 

/KSTG 

/L8CNT 

/LINT0T/MA4 

/  M8  AL 

/NTPU 

/P 

/PC 

/PDC 

/PD r  7pwr  — 

/PHT6 

/pffrr- 

#whilc 

/PHTCHC 

/PP  ~ 

#  PS  I  /PSID 

/PSIE 

/PTIN 

/PXPC 

/PXPP 

/PXPS 

/°YPP  / PYPS 

/PZPC 

/PZPP 

/PZPS 

/P11 

/PI  2 

/  P21  /P  22 

/  P2  3 

/P31 

/P32 

/P33 

/3 

,QC  ,QOC 

/QDS 

/QMIN 

/OP 

/<JS 

/R 

/■PC  /  R  DC 

/  RD  S 

/RET 

/RHO 

/RP 

/RS 

3/SIGMAY/SINPSD 

'/SINTH 

D/SPHI 

/SPSI 

/STHE 

/  S  XGZ  G! 

/TBS  /TDIFF 

/  TDT 

/TEMPI 

/TEMP  2 

/TEMP3 

/THETA 

D/THETAE/THETL 

/TIMC 

/TIMP 

/TIMS 

/TLIM 

/  TST 

X/T$MMIN/TTT 

/  UM  A 

/VIXC 

/VIXP 

/VIXS 

/Virc 

/VI YS  /VIZC 

/VIZP 

/VIZS 

/  VMA 

/VN 

/VPRIX 

/VPRIZ  /VOW 

/  VRXP 

/  V  RY  P 

/  VRZP 

/VS 

/VW 

/VWS  /VX 

/VY 

/VI 

/  W 

/WAC3 

/WAces 

/ WD  AC  3  / W  D A  C  9  S 

l/WDS 

/WS 

/WIN 

/  WMA 

/  WPS 

/  WPN  /WPS 

/Vis 

/XC 

/XDC 

/  X  DD  C 

/XDDS 

/ XD IN  / XDP 

/  XD  S 

/XEL 

/  XET 

/XIN 

/XIXX 

/XIZZ  / X<  N 

/XL  AM 

/XLAMM 

/  XL  CM 

/XLCP 

/  XL1 1 

/ XL  1 3  /XL21 

/  XL22 

/XL23 

/  XL  31 

/  XL3  2 

/XL33 

/ XM A 1 0  /XM ASS 

,  XM  A  X 

/XMAY 

/XMAZ 

/XMAVX 

/XMAVY 

/ XM  C  / XMF  X 

/  XMFY 

/XMFZ 

/  XMSP 

/  XMT  X 

/  XMTY 

/XMXACS/X.MY  ACS 

i/XMZACS/XNLOAF 

:  /  X  N  U 

/XNV1 

/  XN  V2 

/XP  ,XPC 

/  XPDC 

/  X  PD  D 

/XPDS 

/XPP 

/XPS 

/YC  / Y  DC 

/  YDDC 

/YDDS 

/  YDG 

/YDIN 

/  YD° 

/YEL  /YET 

/YMAID 

/YP 

/  YP  C 

/YPDC 

/  YP9D 

/Y P?  /YPS 

/  YS 

/ZC 

/  ZDC 

/ZDDC 

/  ZDDS 

/ZDIN  /ZDP 

/ZDS 

/ZEL 

/ZET 

/ZIN 

/ZMAID 

/ZPC  /ZPOC 

/  IP DD 

/ZPDS 

/ZPP 

/ZP$ 

/ZS 

/  ZZ  2  / Z  Z3  - 

/  ZZ 4 

/  Z  Z5 

/  ZZ5 

/  ZZ7 

/  ZZ8 

/ADE11  / A  DE 1 2 

/ AD  El  3 

/ A  DE  21 

/ ADE22 

/ADE23 

/ ADE31 

/ AD  123  /ALDC 

/ALDCL 

/ ALPHFl 

./ALPHYL 

/  ALP  C 

/ALRC 

/ALPFP  /ALRPY 

/  ALTC 

/ ALTCL 

/AMAX 

/ A®E  ST 

/AYEST 

/APY  / APZ 

/  3  E  0  R 

/BETA 

/3ETASH 

/9ETASW 

/  C  AO 

E/CASF  /CNLIM 

/CUN 

/CUE 

/CUD 

/OELPA 

/ DELR  A 

/OcLPL  /DELOL 

/DELYL 

/Dlam 

/ DL  AMD 

/D  VN 

/  DVE 

/ENSCS  /ESSCS 

/ERALT 

/GACPC 

/GADOPC 

/  G  AD  DPS 

/ GDDDP 

/GADDYC/GACOYS 

l/GDDDY 

/GAMMAD 

'/  GAMMAT 

/  GDD  C 

/ GDDCL 

/GD?C  /GDYC 

/GDPX 

/GDPY 

/GDPZ 

/  G  DTC 

/ GDTCL 

/ H A  3 S  /HOOT 

/  HN  A  V 

/  H  PC  5 

/  IA  SW 

/ 1  ATM 

/  iesw 

/ I F IL  /IGAIN 

/  IL  AG 

/ILAST 

/IR REF 

/IRUP 

/ISCS 

/ 1 S  HD  /ISTER 

/  J  E  S  W 

/N  RE® 

/PHIAL 

/ P  HI CMD 

'/PHIGD 

3/PN  /PE 

/PD 

/POS 

/“UN 

/PUE 

/PUD 

/OXGD  /»KGT 

/  RUX 

/RUY 

/RUZ 

/TEJ 

/TG 

/TlAMD  / TPD 

/  TR  TN 

/UPC 

/URC 

/UYC 

/  VHS 

/ VX  3  /VYS 

/  VZS 

/WAP 

/ X  F  A IL 

/YFAIL 

/ZFAIL 
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TABLE  A-l.  CODE  LISTING  FOR  THE  DATACOKR  SUBROUTINE  (Continued) 


COMMON 

XMT  / XMTD 

/  YMT 

/  YMTD 

/2NT 

/INTO 

END 

COMMON  PACK 

**  * 

COMM ON /OTHERS/ AC 

✓  ACM 

/ACPX 

/ACPY 

/ACPI 

/  ACX 

/ACY 

A 

1  AC  Z 

/APK2/14) 

/APC 

/  AP  FB 

/AYC 

/  AY  FB 

/CGE 

/  CLDR 

A 

2CCSTC 

/DGAMC  /DNVC2/6) 

,EPP 

/EYP 

/FEEP 

/GAMMCS/GAMMG 

A 

3IACFLG 

/I APFLG/IBO 

/ 1 D83 

/ 1 D34 

/IDB5 

/ 1 D  86 

/  IDI3 

/  IDI4 

A 

4IDI5 

/IDI6  /  ICOF 

/ IGA  FLG/ IPOFL 

/P  MX 

/P  MY 

/  PMZ 

/GAP 

/ 

5SINGC 

/SINPC  / SINTC 

/TSSM 

*  TG  SM  S 

/TTIMP 

/ TV  AC  F 

/TVACP 

/Til 

A 

I6T1  2 

/T13  /T21 

,122 

/T23 

/T31 

/T32 

/T33 

/VG 

7VIS 

/VTS  /XMfP 

*  X  MS  AP 

/YIN 

/HANKC20D 

.ircv 

•  ■ 

J 

DINENSKN  DMG(64)/PffG<64) 

DIMENSION  DCGC54)/PCG(54) 

DIMENSION  A$(6/6)/SAS(6/6)/AA(6/6)/AAT(6/6)/ASl(6/6) 
DIMENSION  AL(6/6)/AM(6/6)/AM1(6/6)/A£(6/1)/AN(5/1) 
DIMENSION  DC<6/1)/AP<6/1 > / AR (6/1 ) / AA PC6/1 > 


DIMENSION  X(6/12)/'AOSK<12)  /IHL0<6) 

DATA  DMG/5H 

THE/5H 

PHE/5H 

P  S  E  /  5  H 

XE/5H 

YE/5H 

ZE/ 

3H 

THC/5H 

PHC/5H 

PSC/5H 

XCC/5H 

YCC/5H 

ZCC/ 

5  H 

THDC/5H 

PHDC/5H 

PSDC/5H 

XDCC/5H 

YDCC/5H 

ZDCC/ 

5HTHDDC/5HPHDDC/3HPSDDC/3HXDDCC/5HYDDCC/5HZDDCC/ 

5  H 

AA11/5H 

AA12/5H 

AA1 3/5H 

AA14/5H 

AA15/5H 

A  Al  6  / 

5  H 

AA21/5H 

AA22/5H 

AA22/5H 

A  A  2  4  /  5  H 

AA23/5H 

AA26, 

5  H 

AA31/5H 

A  A  3  2  /  5  H 

AA33/5H 

AA34/5H 

AA35/5H 

AA  36  / 

5H 

AA41 / 5  H 

AA42/5H 

AA43/5H 

A A44 / 5  H 

AA45/5H 

A  A  4  6  / 

5  H 

A  A  5 1  /  5  H 

A  A  5  2  /  5  H 

AA52/5H 

AA54/5H 

AA55/SH 

AA  56/ 

5M 

AA61/5H 

AA62/5H 

AA63/5H 

AA64/5H 

A A  65 / 5  rt 

AAo6  / 

5  H 

IXX/5H 

IYY/5H 

IZZ/5HIMASS/ 

DATA  DCG/5H 

A  S  1  1  r  5  H 

AS12/5H 

AS13/5H 

AS14,5rt 

AS15/5H 

AS1  6/ 

5M 

AS21/5H 

AS22/5H 

AS23/5H 

AS24/5H 

AS25/5H 

A  S  2  6  / 

5h 

AS31/5H 

AS32/5H 

AS33/5H 

AS34/5H 

AS35/3H 

AS36/ 

3M 

AS41 ,5H 

AS42/5H 

•S43/SH 

AS44/5H 

AS45/5H 

AS  46  / 

5  H 

AS51/5H 

AS52/5H 

AS  53/ 5  H 

AS  54/ 5H 

AS55/SH 

A  S  5  6  / 

5  H 

AS61 /5H 

AS62/5H 

AS63/5H 

AS64/5H 

AS65/5H 

AS  66  / 

5H 

0  C 1 1  /  5  -t 

DC21  /5H 

DC  31 / 5  H 

DC41/5H 

DC51/5H 

DC  61  / 

5  H 

AP1 1 /5H 

AR21/5M 

AP31/5H 

AR41/5H 

AR51/5M 

AR61  / 

3  H 

AE11/5H 

AE21/5H 

AE31/3H 

A  E41 / 5  H 

AE51/5H 

AE61/ 

EQUIVALENCE 

(HANK (1 ) 

/ X  PC  DTC  ) 

/ <HANK(2>/YPDDTC)/ <HANK(3)/Z® 

IFCIPRNT.EQ.1)  RETURN 

INITIALIZATION  DONE  ONLY  CNC E- IN TF * V AL 


IFCJJJ.S0.1)  GO  TC  5 
KKK  =  0 

S  CONTINUE 
<<<=<<<♦1 
C 

C  NO.  EX^I^ENTAL  DATA  »C  INTS -I  NP  UT 

£ 

N  =  CIN(?«f«) 
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tart,b  A-l.  CODE  LISTING  FOE  THE  DAXACOER  SUBROUTINE  (Continued) 


WRITE<6/111>  KKK/N 

111  F0RMAT(/55X/SHINTERVAL/3X/I3/lH-/I3//> 
WRITS<6/1112)  TIMt 
1112  F09MAT<56X/'TIME'/G15.7> 

C 

C  COVARIANCE  MATRIX/ INITIAL-INPUT 

C 

IF (KKK.EQ.1)  GO  TO  15 
00  15  1=1/6 
DO  15  J  =  1 / 6 
AS(I/J)=3AS(I/J) 

15  CONTINUE 

I F ( KKK • G  T • 1 )  GO  TO  10 
OC  40  1=1/6 
00  40  J=1/f 
AS (I/J )=0. 

40  CONTINUE 

AS(1/1)=CIN<393) 

ASU/2)  =  CIN(397) 

AS (3/ 3 )  =  C I ^ ( 39  6) 

AS(4/4)=CIN<395) 

AS(5/5)=CIN(394) 

AS<6/6)=CIN(393) 

10  CONTINUE 
J  J  J  =  1 

U70  VALLES  TAKEN  AS  COMPUTED 
C  AND  EX°E  5IMENT  4L  DAT  A-INPUT 

r 

THS=THE7 A 
PMc=PH I 
PS  5  =  PS  I 

7HE=7H?*CIN(392) 

PHE=°HE+CIN(3?1) 

pSE=PSc+CIN(390) 

THC=TH-TA 
PH  C  =  PH I 
PSC=pSI 
THDC=TrtETAD 
pMOC=PH  I  D 
PSDC=pSID 
TH0DC  =  '3DC 
PHODC=pCC 
PSDDC=*CC 
IMAS  3  =  W 
GS7=CIN(  21) 

IXX=XIXX*G3T 
IYY=XIYY*S3T 
IZZ=XIYY*G37 
X5  =  ATC  +  CIN(3i<5) 

Y-=tTC*CIN(309) 
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TABLE  A-l.  CODE  LISTING  FOR  THE  DATACORR  SUBROUTINE  (Continued) 

'TFSffTCKllTCTfT) 

XCC=ATC 

YCC=9TC 

ZCC=GTC 

XDCC=A11*VX+A1 2*VY+A1 3*VI 

YDCC=A21*VX+A22*VY+A23*VZ 

ZDCC=A31»VX+A32*VY+A33*VZ 

X00CC=XP0DTC 

YODCC=YFDDTC 

ZDDCOZP3DTC 

SRORR  MATRIX 


AE(1#1)=THE-THC 
AE(2r1 )=PHE-PHC 
AE(3^1)=PS=-PSC 
A£<4,1 )=XE-XCC 
AE(5,1)=YE-YCC 
A£ (6,1 )=2E-ZCC 

INERTIAL  TO  MISSILE  MATR IX“A3 

CTH  C  =4  OS (THC ) 

5TMC=SIN<THC) 

CPHC=C0S (PHC) 

3PhC=?IKTPHC)  ‘ 

C  PS  C -C  0  5  (PSC) 

SP3C=SIN(P3C) 

AJ11=CPSC*CTHC“S?SC*SPHC*STHC 

A?12-CT,HC*SPSC,fSTHC*S0RC*CPSC 

AP13  =  -STHOCPHC 

AE  21 =-SF3C*C?HC 

»S22=CPHC*CPSC 

AB  Z3  =SPH  C 

».P31=CPSC*STHC  +  CTMC*SPHC*SPSC 
A532=STHC*3PSC-CThC*3PMC*CPSC 
Ab33=CTHC*CPHC 


SENSITIVITY  MATRIX 

A  A (1 ,1  )* 

AA (1 ?2)~ 

AA<1,3>= 

A A (1,4)= 

AA(1,5)= 

AA(1/6)= 

;ac  :,d* 

AA (2/2)= 

AA(2,3)= 

AA<2,4>= 

4A<2,5>= 
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TABLE  A-l.  CODE  LISTING  FOR  THE  DATACORR  SUBROUTINE  (Continued) 

A'A<2*6>  = 

AA ( 3  *  1  )  = 

AA(3 *2>* 

AA(3/3)= 

AA(3*4)= 

AA (3*5  )  = 

AA(3*o)= 

AA (4/1  )  = 

AA<4 *2)  = 

AA(4*3)= 

AA(4*4)= 

AA(4*5)= 

AA<4  *6>  = 

AA(5*1)= 

AA<5*2)= 

AA(5*3)= 

AA(5*4)= 

AA(5*5)= 

AA(5*6)= 

AA ( 6* 1  )  = 

AA (6*2)= 

A  A  (  o  *  3  )  = 

AA(t*4)= 

A  4  (6  *  5 )  = 

AA ( 6*  e  )  = 

IAIN  5PCJP  P 9  I  NT  OPO  Eft 


=  *G(  1)  =  THS 
<>«',(  2 )  =  3H  E 
PI-  3  C  3 )  =  °S  E 
=>K3(  4  )  =  <  E 
°*3(.  5  )  =  YE 
3VG(  6)=3E 
P*G(  7 )  =  THC 
Of'EC  3 )  =  PHC 
P"3(  9  )  =  P  S  C 
®F!G(10)*XCC 
°¥G  ( 1 1  )  =  YCC 
0*3(12)= ZCC 
(1 3)=THJC 
=¥3(14)=PH0C 

°*3 (15)=osoc 

°*j(16)=X0CC 
Pf-G  ( 1  7)  =  Y0CC 
°*3(1S)=ZDCC 
°*G(19)='rH0CC 
°*G(ZO)==HDOC 
•>*3(21  >=?S3DC 
o*3( 22)=XD7CC 
3¥3(23)=Y03CC 
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TABLE  A-l.  CODE  LISTING  FOE  THE  DATACORR  SUBROUTINE  (Continued) 

f»«3C24>-ZDDCC 

MNG(25)=4A<1,1) 

®M3<26)=  AA  (1*2) 

®MG(27)=AA(1/3) 

PM  3(28)  =  AAd,4) 

PM3(29)=AA(1*5) 

PMG ( 30)  =  AA (1 *6) 

-PMG(31)=AA(2*1) 

PMG(32)=AA(2*2) 

®PG (33 ) =AA(2*3) 

PM3(34)=AA(2/4) 

PMG (35 )  =  A A (2*  5) 

®MG(36)=AA(2*6) 

PMG ( 37)=  AA  < 3  /I  ) 

PMG (33)=4A(3*2) 

PM3(39)  =  AA  (3*3) 

PMG (  40 )  =  A  A (3*4) 

PMG(41 )  =  AA  (3*5) 

PM3(42)=AA(3*6) 

?mg(43)=aa<4,1) 

PMG(44)=AA(4*2) 

?MG(A5)=AA (4*3) 

?MS(4£)=AA(4*4) 

°M8(47)=AA(4*5) 

°M3(46)=AA(4*6) 

?M3(49)=A4(5*1) 

?MG(50)=AA(5*2) 

®M 8 ( 51 )  =  AA  (5*3) 

?M3(52)=AA(5*4) 

°MG (53)=AA(5*5) 

®MG(54)= AA(5*6) 

°M3<55)=AA(o*1) 

PMG(56)  =  AA  (6*2) 

PMG(57)=AA(6*3) 

°MG  (  5  8  )  =  A  A  (6  *  4  ) 

PM8(59)= AA (6*5) 
p»3(^C')  =  AA  (6/6) 

PM3Co1)=IXX 
®MG  (  ?2 ) = I Y Y 
PM3<  33)  =  III 
=>"G(64)  =  IMASS 

25  FC?  MAT  (OH  *  6  (  1  X  /  A  5  *  31  c  .  7  )  )  ) 

WSIT=(6,2Z5)  (DM3(T)*pyG(I)*I=1*64) 


CCMPUT3  AA  TpA\Sp03E=A AT 

IKi:0N(7G).£0.C)  30  TO  30 

N  =  6 

■'•M 1  =  N“  1 
OC  20  T=1*NM1 
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TABLE  A— 1.  CODE  LISTING  FOR  THE  DATACORR  SUBROUTINE  (Continued) 


IP1=I+1 

00  2 U  J=IP1,N 

TM?1=A 4(1/ J) 

AA(I,J)=AAT{J,I) 

20  AAT ( J , J) =TMP1 

wRITE(o,1002)  A  A  T  <  3 , 1 ) 

1002  PC  RMAT ( 1 X, 31  5.7) 

30  CONTINUE 

If (K<K.tD.1>  GO  TO  35 
C 

C  COMPUTE  INVERSE  AS  =  A  S 1 

C 

DO  50  I=1/i 
DO  50  J  =  1,6 
XCI,J)  -  AS<I,J) 

50  CONTINUE 
MN1=6 
N=  6 
N9=N 
MS  =1 
IC=1 
I C  =0 
IS  =  1 

call  SEaOMI(X,N,N?,HS,MN1,D,R,5,WORK,IHLD,IC,ID,IS) 

DO  6  0  j  =  1,o 
DO  6 J  <=1,6 
AS1(J,<)=X(J,K) 

6 J  CONTINUE 
15  CONTINUE 

IFCKXK.OT.I)  GO  TO  5 15 
JO  5  5  I  =  1,«> 

DO  55  J  =  1 , 5 
AS1 <I,J)=0. 

AS 1 <1,13=1  . / C IN ( 393) 

AS  1  (2, 2)  =1  .  / Cl N( ?°7) 

A  S 1 ( 3, 3 ) =1  ./CIN(  396) 

A  S 1 (4, ) =1 ./CIN<395) 

A  SI (  5 , 5  )  =  1  ,/CIN(3«n) 

AS1 (0,5) =1 ./Cl N (393) 

55  CONTINUE 

APIT5(14,40C) 

400  f OR»AT (1 x,  •  55* ) 

r 

'C  SrNSITIV:TY-T*CCVA^IANC-1  *ATPIX=AL 

C 

M  =  5 
N=  6 
3  =  6 

DC  70  1=1, M 
DO  70  J  =  1,=> 

4L (I, J  )=  ). 
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TABLE  A— 1.  CODE  LISTING  FOR  THE  DATACORR  SUBROUTINE  (Continued) 
0  0  70  K*1/N 

70  ALU,J)  =  AL<I,J)  ♦  AAOC,I)  *  ASKK/J) 

WRIT  6(14*401 ) 

4u1  FORMAT (1 X, ' 70 ' ) 

AL*SENSITI VITY  MATSIX=AM 

n-6 

N=o 
P=  6 

DO  80  1=1, M 
DO  80  J  =  1,P 
AM(I,J>=0. 

DO  80  K=1  ,N 

SO  AM<I,J)=AM(I,J)  ♦  A  L  ( I ,  K  }  *  A  A  (K  ,  J  ) 

WRITE(14,402) 

-*02  FORMAT  (1  X,  '  80'  ) 

AL*  A  E  MA  TR IX  =  AN 
N  =  6 
?=1 

DO  90  1=1, M 
DO  RO  J  =  1,=» 

AN(I,J)=0. 

DO  90  <=1,N 

90  AN(I,J)  =  AN(I,J)  +  AL<I,K)  *  AE(K,J) 

WRITS(14,4G3) 

4C3  e0RMAT(1X,,90*) 

COMPUTE  INVERSE  AM  =AM1 

DC  100  1=1,6 
DC  100  J  =1 ,6 
X(I,J)=AM(I,J) 
luJ  CCNITINJe 

WRIT  EC  14,404) 

404  F0RMATC1X,  MOO 
MN1  =  6 
\=0 
NO=N 
MS  =  1 
I  C  =  1 
I D  =  C 
I  S  =1 

CALL  S  =  S3MICX,N,\8,MS,MM,D,R,E,W0RK,IrtLD,IC,ID,IS> 
DO  110  J  =1 , e 
DC  110  K=1,t 
110  A  M 1 ( J , < )  =  X  (  J  ,  A  ) 

«.p:TEC14,405  ) 

•♦‘•E  FC3MAT(1X,'1l6‘r) 
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TABLE  A-l .  CODE  LISTING  FOR  THE  DATACOKR  SUBROUTINE  (Continued) 

C  AMI  *  AN  P AT9IX=DC 

C 

<1  =  6 
N=6 
P=1 

DO  120  1=1, M 
DC  120  J=1,P 
0  C  <  I ,  J  )  =  0  . 

DO  120  K  =1  ,  N 

120  DC  (I,J)=QC  (I,J)  ♦  Af*1  (I/O  *  AN(K,J> 

C 

c  s£msitiv:ty*dc=ap 

c 

vi  =  o 
N  =  M 

D  =  1 

DO  130  1=1, M 
DO  130  J=1,P 
AP(1,J>0. 

DO  1  30  K  =1 , N 

120  AP(I,J)  =  AP(I,J>  +  AA(I,K)  *  DC  (  K  ,  J  3 

s, 

:  Ac”  AP  1»TPIX=C 

C 

N?=6 

SC=1 

dc  no  j  =  i , n c 

DC  no  2=1, N? 

140  4A?  <I,J)  =  A€(I,J>  -  4P ( I , J ) 

N=CI\<3W) 


0*3  ”T  =  A  9 

AS(1,1)*AA'<(1,1)*AAR(1,1)/N 
A3(1,2)=AAP(1,1)*AAR(2,1)/3! 
A3(1,3)=AAS(1,1)*»ARC2,1)/N 
Ao(1,4)=aas(1,1)*AAR(4,1)/N 
•S«1,?)=AA?f1,1)*RAR(5,1)/N 
AS(1,6)=AACC1,1) *AAR(t,1)/N 
4  S  <  2 ,1 )  =  AA5(2,1)*AAR  (1 ,1  )/N 
4S(2,2)=»AA(2,1)*AAPC2,1)/\ 
AS(2,7)=AAR(?,1)*AAR(3,1)/N 
AS(2,4)=AA1(2,1)*AAsC4,1>/M 
AS(2,?)  =  AAA(2,1)*’AAP(5,1)/N 
AS(2,6)=AA3(Z,1)*AAS(6,1)/N 
aS<3,1)=AAP(3,1)*AAP(1,1)/N 
AS(3,:)  =  A4A(?,1)*AAP<Z,1)/'< 
4j(3,7)=*A‘><3,1)*A4A(:,1)/M 
A5(:,h)=AAR(3,1)*A4p(4,1)/4 
AS(I,5)“AA3(3,1)*AAP(5,1)/S 
4«(3,6)=»A5(;,i>*aap(6,1)/N 
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TABLE  A— 1.  CODE  LISTING  FOR  THE  DATACORR  SUBROUTINE  (Continued) 

AS<4/1)*AAfl(4/1)*AAR(1/1)/N 

AS(4/2)=AAR(4/1)*AAR(2/1)/N 

AS(4,3)=AAR(4/1)*AAR(3/1)/N 

AS(4/4)=AAR(4/1)»AAR(4/1)/N _  _ 

AS(4/5)=AAR(4/1)*A4R(5/1)/N 
AS(4/6)=AAR(4/1)*AAR(6/1)/N 
AS(5/1)  =  AA:RC5/1)*AAR(1/1)/N 
AS(5/2)=AAR(5/1)*AAR(2/1)/N 
AS(5/3)sAAR(5/1)*AAR(3/1)/N 
AS(5/4)=AAR(5/1)*A4R(4/1)/N 
AS(5/5)=AAR(5/1)*AAR(5/1)/N 
AS(5/6)=AAR(5/1)*AAR(6/1)/N 
4S(6/1)=AAR(6/1)*AAR(1/1)/N 
AS(6/2)=AAR(6/1)*AAR(2/1)/N 
AS(6/3)=4AR(6/1)*AAR(3/1)/N 
4S(6/4)  =  AAR(;>/1)  *AAR(4/1)/N 
4S(o/5)=AAR(6/1)*AARC5/t)/N 
AS (6/6)=AAa(6/1) *  A  AR (6/1 ) /  N 
)  C0R,VAT(S7X/17HC0NVER5ENCE  GROUP) 

RESIDUAL  MATRIX-ONE  ITERATION  AHEAD  ERRCR  MATRIX 


AR (1  ,1 > 
AP(2/1  ) 
A  R  (  3  / 1 ) 
A  P,  (  4  / 1  ) 
AS(5/1) 

A  R  (  ft  / 1  ) 


4  A  3  (  1  /  1  ) 
A  A  R  (  2/  1  ) 
4A4C3/1) 
A  A  R ( 4/ 1 ) 
AA3<  3/1  ) 
AA3C6/1) 


CONVERGENCE  GROUPS  PRINT  ORDER 


P  C  G  (  1) 
aCi(  2) 
RCGC  3) 
aCG(  4) 

°  C  3 (  5) 
3  C  3  (  fc) 
°C3(  7) 
aC  3  (  i) 
aC'j(  9) 
3  C  G  ( 1  0  ) 
PC  5 ( 1 1  ) 
°  C  G  <  1  2  ) 
°C  3(13) 
PC  G ( 1 4  ) 
PC3C1S) 
aC  G ( 1 6  ) 
PC  G (17) 
3  C  3  ( 1  ) 
3  C  3  (  1  9  ) 
°C  i(2C) 


*AS(1/1) 
=*5(1/3) 
=  A  S  ( 1  /  3 ) 
=AS (1/4) 
=A$(1/E) 
=AS (1/t) 
=Ab (3/1 ) 
=40(2/2) 
=AS(2/3) 
=  A  S (2/4) 
= AS(2/5) 
=  AS (2 /C ) 
=AS (3/1 > 
=  AS (3 / 2 ) 
=AS (3/3) 
=AS(3/4) 
=  A  3  ( 3  /  5  ) 
a  AS (3/6) 
*AS(4/1) 
-AG (4/2) 


TABLE  A-l.  CODE  LISTING  FOE  THE  DAIACOBE  SUBROUTINE  (Continued) 


PCG(21)=AS(4,3) 

BC3(22)=AS(4,4) 

PCG<23)*AS(4,5) 

PCG(24)=AS(4,6) 

PC3(25)=AS(5,1) 

3C3(26)=AS (5,2) 

°C  3 ( 27)  =  4S (5,3 ) 
eC3(2*!)=AS<5,4) 

PC3(29)=AS (5,5) 

°C  3( 30) =AS (5,6) 

PCG(31)=AS(6,1) 

°C  G ( 32 )=  AS ( 6, 2) 

°C3(33)=4S (6,3) 

PCGC34)= AS(6,4) 

PCG(35)=AS(6,5) 

PC3(3e>)=  AS  (6,6) 

=>C3(37)  =  DC  (1,1 ) 

°C3( 32  )  =  D  C  (2,1 ) 

PC3(3<?)  =  OC(2,1) 

®C3(4G)=DC(4,1> 

PC  3(41 )  =  DC  (5,1 ) 

?CG(42)  =  DC  (6,1) 

3C3(43)=AR(1 ,1  ) 

°C  G( 44 )  =  A R  (2,1 ) 

?C3(45)=AR(3,1 ) 

3CG(46)=AR(4,1) 

°C3<47)= AP(5,1  ) 
dC3(4?)=A3(6,1) 

PC3(49)=A£(1,1) 

PC3(50)=A'(2,1) 

°C3(51)  =  A*  (3,1) 

°C3(52)=A£(4,1) 

3C3(53)=At(5,1) 

°CG<54)=  AE  (6,1 ) 

C 

C  WRITS  CCVARIANC,DC,R£SIOUAl,  =  RRC.R 

r 

s* 

WRITE (a, 500) 

WRITE(  6,  225)  (DC3(I),PCG (I),  1  =  1,54) 
DC  150  1=1,6 
DO  150  J  =1 ,6 
15  )  3AS(I,J)  =  A S  (  I , J  ) 

RETURN 

END 
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APPENDIX  B 

MODIFICATIONS  TO  BOUT  (Boost  Output) 
AND  RVOUT  (Re-entry  Vehicle  Output) 
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APPENDIX  B 


MODIFICATIONS  TO  BOUT  (Boost  Output) 

AND  RVOUT  (Re-entry  Vehicle  Output) 

The  changes  to  the  BOUT  and  RVOUT  subroutines  of  the  U70  simulator 
needed  for  the  Tactical  Ballistic  Missile  Trajectory  Study  are  listed  in 
Table  B-l. 


TABLE  B-l.  MODIFICATIONS  TO  BOUT  AND  RVOUT 

PRM(70)  changed  to  PRM  (72) 

RPG  (162)  added 

DIMENSION  PR(30) ,PRE( 18) ,  PRM(72) ,PRG(78) ,RPG(162) 

Both  DO  loops  changed  from  '1,70'  to  '1,72' 

725  DO  730  IMX  -  1,72 

290  CONTINUE 

IF(IMISS.EQ.O)GO  TO  310 

Added  after  PRM(70)  -  WACB 

PRM( 71)  -  0 

PRM(72)  -  0 

551  CONTINUE 

PRG(22)-  PRG(22)*RTOD 
PRG(23)»  PRG(23)*RTOD 
PRG(28)-  PRG(28)*RTOD 
PRG(29)-  PRG(29)*RTOD 
PRG(34)-  PRG( 34 ) *RT0D 
PRG(35)-  PRG(35)*RTOD 
PRG(41)-  PRG(41)*RT0D 
PRG(64)-  PRG(64)*RTOD 
PRG(65)-  PRG(65)*RTOD 
PRG(75)»  PRG(75)*RTOD 

RPT(I)  -  PRG(I)  added  after  PRG(75)  -  PRG(75)*RTOD 

RPG(4)-PRG(22) 

RPG(52)-PRG(23) 

RPG(10)-PRG(28) 

RPG(58)-PRG(29) 

RPG(16)-PRG(34) 

RPG(64)-PRG(35) 

RPG( 137) “PRG( 4 1 ) 

RPG(130)-PRG(64) 

RPG(151)-PRG(65) 

RPG(129)-PRG(75) 
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TABLE  B-l.  MODIFICATIONS  TO  BOUT  AND  RVOUT  (Continued) 


WRITE  (14,26)  statements  added  before  the  WRITE  (6,2)  statements 

IF(ICON(52).EQ.O)  GO  TO  515 
WRITE( 14,25)  PTIM 

515  CONTINUE 

25  FORMAK615.7) 

C  WRITE(6,2)  PTIM.TSTT , DT , IRJU , IFAIL 

IF(ICGN(52).EQ.O)  GO  TO  516 
WRITE( 14,25)  TIMC 

516  CONTINUE 

C  WRITE(6,2)  TIMC, TST,DT, IRJU, IFAIL 

IF(ICON(52).EQ.O)  GO  TO  533 
WRITE( 14,26)  (PR(I), 1-1,30) 

533  CONTINUE 

WRITE(6 ,5)  (DR(I),PR(I), 1-1,30) 

IF(ICON(52).EQ.O)  GO  TO  534 
WRITE(14,26)  (PRE( I), 1-1,18) 

515  CONTINUE 

26  FORMAT(6(E15.7) ) 

WRITE(6 ,2)  DRE(I),PRE(I), 1-1,18) 

IF(ICON(52).EQ.O)  GO  TO  546 
WRITE( 14,26)  (PRM(I) ,1-72) 

546  CONTINUE 

WRITE(6 ,5)  (DRM( I), PRM(I), 1-1,70) 

IF( ICON(52) .EQ.O)  GO  TO  553 
WRITE(14,26)  (RPG(I), 1-1,162) 

515  CONTINUE 

WRITE(6 ,5)  (DRG( I) ,PRG( I) ,1-1 , Ml) 
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TABLE  B-l.  MODIFICATIONS  TO  BOUT  AND  RVOUT  (Continued) 

PRG<78)  tff~  - 

IF(IC0N(52).EQ.G)  SO  TO  220 
RR6<  1>  »  PRG(t) 

RRG<  2)  a  prg<2) 

RRS<  3)  a  PRG (3) 

RRG(  4)  a  prg(22) 

RfttH  5Y  ~  PftG(^) - 

RRGC  6)  a  PRG(6) 

RRGt  7)  a  PRffC?) 

RRG (  S)  a  prg(3) 

RRG<  9)  *  PRG (9) 

RRG(IO)  a  PRG<28) 

RRG<11)  a  prg<11) 

RRG<12)  a  PRG (1 2 ) 

RRG<1 3)  *  PRGCTir 
RRGC14)  a  PRG(1 4) 

RRG<15)  a  PRG<T5) 

RRG ( 1 6)  a  prg(34) 

RRG(17>  a  PRGdTF  - 

RRG <1 8)  a  PRG (1 8 ) 

RRG(19>  a  0 
RRG  <  20)  =  0 
RRG  <21 )  =  0 
RRG(22)  *  o 

RRS<23)  a  0 -  -  -  -  - 

RRG  <  24)  a  0 
RR3(25)  a  Q 
RRG(26)  a  0 
RRG(27)  a  0 
RRG (28)  a  0 

RRG( 293  s'PirSttO) -  -  ' 

RRG(30)  a  0 
RRG<31 )  =  0 
RRGC32)  »  0 
RRG <  33)  a  a 
RRG<34)  =  0 

RRG<35>  a  Q- - — 

RRG ( 36)  *  o 
RRG<  37)  =  0 
RRG( 33 )  a  o 
RRG<39)  a  0 
RRG(40)  *  o 

RRG  <41 )  a  - - - 

RRGC42)  »  o 
RRG (43 )  *  0 
RRG(44)  a  o 
RRG(45)  a  0 

INOGC  *  I APFL6  *  10*IGAFLG  ♦  100*IACFLG  ♦  1GOO*ISE<t 
RINOSC  »  INDSC  ♦  99* J>5  ♦  0.000005 
RRG ( 46 )  a  PRG<76 ) 
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TABLE  B-l.  MODIFICATIONS  TO  BOUT  AND  RVOUT  (Continued) 


RRG(47> 
RRG(48) 
RRG<49) 
RRG(SO) 
RRG*51> 
RRG(52) 
RRG(53) 
RRG ( 54) 
RRGC55) 
RRG ( 56) 
RRS<57)“ 
RRG(58) 
RRS(59) 
RRG ( 60) 
RRG(61) 
RRG (62) 
RR€(63> 
RRGC64) 
RRG(65) 
IFTEMP 
1  ♦ 


0 

0 

0 

0 

-o - 

PRG(23) 

0 

0 

0 

0 

0 - 


PRGC29) 

PRG<53) 

0 
0 
0 

*  P  R  G  (  3  5  ) 

=  PRG (54) 

=  IALTU*10*(IPUDFL+1)+100*(IUPDFL+1 >+1000*jaSW 
10CC0*(IRCHFG*1) 


RRG ( 66) 


RRG<  67~)  =~ 
RRG(68)  * 
RRG(69)  * 
RRG (70)  = 
RRG(71 )  = 
RRG(72)  = 
RRG(73)  = 
RRG(74)  = 
RRG(75)  = 
RRG(76)  = 
RR3(77)  = 
RRG(78)  = 
TTRG(79> 
RRG( 30) 
RRGC81J 
RRG( 32) 
RRG( 33) 

RRG (84 ) 
RffGT857 
RRG (36) 
RRG(87) 
RRG(88) 
RRG(39) 

RRG ( 90) 
RKG(9t) 
RRG(92) 
RRG(93) 
RRG(94) 
RRG(95) 


TABLE  B-l.  MODIFICATIONS  TO  BOUT  AND  RVOUT  (Continued) 


RRG<96)  =  0 
RRG(97)  *  (T 
RRG ( 98)  =  0 
RRGC99)  *  PRGC69) 
RRG<100)=  0 
— RRGC101)*  0 
RRG<102)*  0 

-RRGT103)*  0 - 

RRG < 1 04)  *  PRG  <4) 
RRG<105>=  0 
RRG<106)*  0 

-  RRG<107>=  0 - 

RRG<108)=  0 

- RRG<109)=  0 - 

RRGd  1 0)  *  PRG<10) 
RRGdll)*  0 
RRG<112)=  0 
RRG (11 3) =  0 
RRG<114)=  0 

-  RRG(tt5)^0 - 

RRG(1 1 6) *  PRG(1 6 ) 

RR3<117)=  0  - 

RRG <11 8)=  0 
RRG (119) -  Q 
RRG(1 20) *  0 

- RRG  (1 21  )  =  PRG  <20 

RRGd  22)  *  PRGC26) 
RRG < 123)*  PRSC32) 


RRG(1 24) * 

PR  G  <47 ) 

RRG(125) * 

PRG<4€) 

RRG< 1 26) * 

PRG<77) 

RRG  (1277*" 

PRG (59) 

RRG  < 1 28) * 

PR  G  <  60  ) 

RRGd  29)  * 

PRG  (75)- 

RRG<1 30) * 

PR  G  <  64 ) 

RRGdrrr* 

PR^(24 ) 

RRGd  32)  * 

PRG  <30 

—  RRG  <13  3)*" 

PRS<36) 

RRG  <1 34)* 

PRG<37) 

RRGU35)=PRS<38) - 


RRG<136>*  PRG<39) 

RR3<137)=  PE5(4T) - 

RRSd38>  =  PRG<42) 

RRG(1 39) *  PRG<44) 
RRGC140) *  PR S <45 ) 

- RRGC  >4  >)*  PRG (49) - 

RRG(142)*  PRG<43) 

— — RRG(T43)*^  PRGT50) 

RRG(1 44) s  PRG<51 ) 

- RR  G  <1 4  3  )g  PR3(32) - 

RRG<1 46) *  PRG<55) 


RRG<148)= 
RRG<149)= 
RRG(1 50) J 
RRS<131^ 


PR  G  <  58 ) 
PRS(61) 
PRG  <62) 
PR  3 (65) 


RRGd  52)  *  PRG  (66) 


RRG  <1 56) 

Rpsdsrr 

RRG<158) 

RR3<160): 


PR6<71 ) 
PR  5  (7?)' 
PRG<73) 
PRG<  7 4) 
PRG<46) 


RRG<1 62) =  0 


PRG<22)*  PRG (c2) *RTOD 

- P*3<73  )*  PRG  <  23h*'RfOO^ 

PRG (28)*  PRG  <28) *RTOD 
PRG (29)*  PRG<29)*RT0D 
PRG(34)=  PRG (34) *RTOD 
PRG<35>-  PRG<35)*RTOD 
PRG<  41 )=  PRG  <41 ) *RTCD 

- PR5<64  >=PRS(u*)*RT0t>- 

PRG<65)=  PRG (65) *RTQD 
PRG (75)  *  PRG<75)*RT0D 
RRG<4)  *  PRG (22) 
RR3*52)-  PRG (23) 
RRG<10)=  PRG (28) 

PRG  <  5  8  >  =  PRG ( 29) 

RRGd  6)=  PRG  (34) 
RRG<64)=  PRG  (35) 
RRG(137)*PRG<41) 
RRG(130)*PRG(64) 
RRG(151)*PRG<65) 

RRGd  29)*  PRG <75 ) 


TABLE  B-l.  MODIFICATIONS  TO  BOUT  AND  RVOUT  (Continued) 


IF(ICON(52).EQ.O)  GO  TO  515 
WRITE( 14,25)  PTIM 

515  CONTINUE 

25  FORMAT(615.7) 

C  WRITE(6 ,2)  PTIM,TSTT,DT, IRJU, IFAIL 

IF( ICON(52) .EQ.O)  GO  TO  516 
WRITE(14,25)  TIMC 

516  CONTINUE 

C  WRITE(6,2)  TIMC, TST.DT, IRJU, IFAIL 

IF(ICON(52).EQ.O)  GO  TO  533 
WRITE( 14,26)  (PR( I), 1-1,30) 

533  CONTINUE 

26  FORMAT(6,(615.7) 

WRITE(6 ,5)  (BT(I) ,PT( I) ,1-1 ,30) 

IF(ICON(52).EQ.O)  GO  TO  534 
WRITE( 14,26)  (PV(I), 1-1,18) 

534  CONTINUE 

WRITE(6,5)  (DV(I) ,PV( I) ,1-1 ,18) 

IF( ICON( 15) .GT.O)  GO  TO  544 
LINC  -  LINC  +  14 
WRITE(6 ,4) 

IF( ICON(52) .EQ.O)  GO  TO  541 
WRITE(14,26)  (RM(I), 1-1,72) 

541  CONTINUE 

IF( ICON(52) .EQ.O)  GO  TO  554 
WRITE( 14,26)  (RG(I), 1-1,162) 

554  CONTINUE 

WRITEC6.5)  (DG(I) ,PG( I) ,1-1 ,NG) 

DIMENSION  PT(30) ,PV(18) ,PM(54) ,PG(120) ,RG(162) ,RM(72) 


RG( 162) ,RM(72)  added 


TABLE  B-l.  MODIFICATIONS  TO  BOUT  AND  RVOUT  (Continued) 


PM(  1) 
PM(  4) 
PM(  6) 
PMt 
PM(1 1 ) 
PM(16) 
PM(1 8) 
PH(1 9) 
PM (20) 
PM<22) 
PM (23) 
PM<24) 
PM(26) 
PM  (27) 
PM(28) 


PM(  1 ) *  GST 
PM(  4)*6ST 
PM(  6)* SST 
PM(  9 7» GST 
PM(11)*GST 
PM(16)*G$T 
PM(18)*GST 
PM(1 9) *G$T 
PN(20)*GST 
PM(22)*GST 
PM(23)*6ST 
PM(24)*GST 
PM(26)*G$T 
PM(27)*GST 
PM(28)*RT0D 


RM(  T) 

RM(  4) 

imc  TT 

RM(IO) 

RW(T33” 

RM (1 9) 

RMC23) 

RH (24) 

RMC2ST 

RM(29) 

RMC30) 

RM (31 ) 

RM(35) 

RM (36) 

RMt37) 


PMC  1) 
PM(  4) 
PMT  67  - 
PM(  9) 

mnr 

PM ( 1 6) 
PM  (18) 
PM ( 1 9) 
PMT  2  01“ 
PM( 22) 
PWC73) 
PM( 24) 
PM( 26) 
PM( 27) 
pm  cm 


PM (29) 

=  PM(  29  r*RT0D - 

RM (3  3) 

3 

PM ( 29) 

PM(30) 

*  PM(30)*RT0D 

RM (39) 

3 

PM ( 30) 

PM (31) 

a  PM(31)*RT0D 

RM (40) 

= 

PM( 31 ) 

PM (32) 

=  PN(32)*G$T 

RMC4T) 

3 

PM(32) 

PM (33) 

a  PM(33) *RT0D 

RM(43) 

3 

PM(33) 

PM(34) 

a  ?M(34)*RT0D 

"RMT44ra 

PM (34) 

PM  (3  5) 

*  PMC35)*RTW~ 

RM(45) 

3 

PM(35) 

PM(36) 

=  PM(36)*RT0D 

RMT46) 

"3” 

PM(36T 

PM(37) 

=  PM(37) *GST 

RM(47) 

= 

PM(37) 

PM (33) 

a  PM ( 38) *RT0D 

RMC497 

3 

PMC38T  ~ 

PM (39) 

a  PM(39)*RT0D 

RM (50) 

3 

PM(39) 

PM (40) 

=  PM(40)*RT0D 

- RMt51 ) 

PMC40) 

PMT4TI 

=  PMC4TT*FTao  ' 

RM(52) 

x 

PM<  41 ) 

PM (42) 

=  PM(42) *GST 

PM (53) 

3 

PM ( 42) 

PM(46) 

a  PM(46) *RTOD 

RM(53) 

= 

PM(46) 

PM(52) 

=  PM( 52) *  RTOD 

RM (67) 

3 

PM( 52) 

TABLE  B-l.  MODIFICATIONS  TO  B(KJT  AND  RVOUT  (Continued) 

PS(120)*DVD-ZNTD 
RGC  t)  *  PG(1>  - 
R6(  2)  »  PG (2) 

RGC  3)  =  PG(3> 

RGC  4)  =  PG(4) 


RG  (5  y 
RGC  6)  = 
RGC  7)  * 
RGC  8)  * 
RGC  9 r  = 
RGCIO)  - 
RGC'H')'  a~ 

RGC12)  » 
RGC13)  * 
RGC14)  * 
RGCt5)  = 
RG  Cl  6)  a 
RGCTTT^a 
RGC13)  = 
RGC197  a 
RGC20)  = 
RGT21)  = 
RGC22)  = 
RGC23)  a 
RGC24)  a 
RGC25)  * 
RGC26)  * 
RGC27)  = 
RGC23)  * 
RGC 2?)= 


a  PG C6) 
a  PG<  7) 
a  PGC8) 
a  PGC9) 
a  PGC10) 
a-pGCttr- 
a  PGC12) 
a  PGC13r 
a  PGC14) 
a  PGC15J 
a  PS  C 1 6) 
PGCT7) 
PGC18) 
PGCT9J 
PGC20) 
PGCztr 
PGC22) 


RGC30) 

RGT31) 

RGC32) 

RGC337 

R6C34) 

RGC35) 

RGC36) 

RGC37) 

RGC3S) 

RGC3?) 

RG (40) 

RGC41) 

RGC42) 

RGC43) 

RGC44) 

RGC45) 

XNDGC  = 

RINOGC 

RG (46) 

RG (47) 

RGC48) 


'  PGC24) 

•  PGt  25) 
PG  C  26) 
PGC28J 
PGC23) 
PGC29) 
PG  C  30) 
PGC3T) 
PGC32) 
PGC33) 
PGC34) 
PG  C  35) 
PG  C  36) 

>  PGC373 
PG  C  38) 
PG  C  39) 
PGC40) 
PGC4t)  ' 
PG  C  42) 
PG  C  43) 
PGC44) 
PGC45) 
IAPFLG  ♦ 

'  INDGC  ♦ 
■  PGC46) 

»  PGC47) 

'  PGC48) 


10*IGAFLG  ♦  1 0C  *1  AC  FLG  +  1000*ISEQ 
99. D5  ♦  0.000005 
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TABLE  B-l.  MODIFICATIONS  TO  BOUT  AND  RVOUT  (Continued) 

RG(49)  *  PG ( 49) 

RG (SO)  =  PG(50) 

RG(51)  *  PG(51> 

R6(52)  *  PS ( 52) 

RG<53)  =  PG<53) 

RG (54)  *  PG(54) 

R6(55)  =  PG ( 55) 

RG(56)  *  PG( 56) 

RG(57>  -  "Pff<  57") - - - - 

RG (53)  »  PG(58) 

RG(59)  =  PG ( 59) 

RG (60)  »  PG( 60) 

RG(61)  »  PG( 61 ) 

RG (62)  =  P G ( 62 ) 

RG (63)  *  PG<63) - -  - 

RG (64)  *  PG(64) 

RG (65)  *  PG(65) 

IFTEHP  =  IALTU  ♦  1  0*  ( I  Pll  D  FL+ 1 )  ♦  1 00* ( IUPD FL+1 >  ♦  1000 
1  J8SW  *■  1 00C0*( I RCHFG+1 ) 

RG (66)  *  P6C66) 

RG (67)  *  PG<67>  -  - 

RG (68 )  -  PG( 68) 

RG (69)  *  PG (69) 

RG(70)  =  PG(70) 

RG (71 )  =  PG(71) 

RG (72)  =  P3(72) 

RG  (  7  3  )  *  PG  (  73) - -  "  "  —  -  "  -  - 

RG (74)  *  PG ( 74) 

RG ( 75 )  =  PG ( 75) 

RG (76)  =  PG(76) 

RG (77)  *  PG(77) 

RG (78)  =  PG (78) 

-RG  -  PG  (  79)" - -  -  - 


RG (80) 

= 

PG ( 30) 

RG (81 ) 

z 

PG(8T) 

RG(82) 

= 

PG( 32) 

RG<83) 

s 

PGC83) 

RG (84) 

= 

P  G ( 84 ) 

- RG(85> 

z 

PG  (  35) 

RG ( 3  6) 

3 

PG ( 86) 

RG(87) 

s 

PG( 87) 

RG (83) 

3 

PG ( 38 ) 

RG (89) 

Z 

PS<39) 

RG (90) 

= 

PG(90) 

-RG(9t) 

— 3* 

psm) 

RG (92) 

3 

PG(92) 

RG(93) 

3 

P3(93) 

RG (94) 

3 

PG ( 94 ) 

RG(95) 

3 

P3(95) 

RG(96) 

3 

PG ( 96 ) 

RG (97) 

S 

PG ( 97) 
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TABLE  B-l.  MODIFICATIONS  TO  BOUT  AND  RVOUT  (Continued) 


RG (93)  = 

PGC98) 

RG (99)  * 

PG<99) 

RG (1 00) s 

PG(IOO) 

RG (101)= 

PGC101) 

R6(102)= 

PG  <1 02) 

R  G  ( 1 Q  5  )  = 

PGClQJ) 

RG (1 04)= 

PGC104) 

RG (105)= 

PGCT05) 

RG (1 06)= 

PG( 1 06) 

R6C1Q7)  = 

PGC107) 

RG (1 08) = 

PG (1 08) 

RG  tT09r)  = 

FG(T09) -  ~ 

R6<110)= 

P3C110) 

RG(111)= 

PG  < 11 1 ) 

RG<1 1 2)= 

PG(11 2) 

RG (1 1 J)= 

PGC113) 

RG (1 1 4)= 

PGC114) 

RGCT15)= 

P5CTT5) - 

RG (1 1 6 )= 

PGC116) 

RG  Cl  17)  = 

PSC117) 

RG (1 1 8) = 

PGC11 8) 

RG (1 1 9) = 

PG  C 1 19) 

RG (120)= 

PG (1 20) 

RG  CT21 )  = 

o  - 

RG<122)= 

Q 

RGC123)* 

IT 

RG<124)= 

0 

RG  Cl  25)  = 

0 

RG (1 26)= 

0 

RGC127)= 

o - 

RG(1 23)* 

0 

RG  129)= 

0 

RG (1 30)= 

0 

R6CT31 ) = 

0 

RG (1 32)= 

0 

0 

RG<134)= 

0 

RGCT35)* 

0 

RG<136)= 

0 

RG<t37)  = 

0 -  -  - 

RG<1 38>= 

0 

RGCT39)® 

0 

RG(14Q)= 

0 

Rorrrrr* 

0“ 

RG (1 42 )= 

0 

JCGCU30* 

0 -  ~ 

R6 (1 44) = 

0 

RG<146)= 

0 

R€C147)= 

0 

RG (1 48) = 

0 

RGC149)» 
RG  <1  50)  = 


R6C152)* 
RGCT5J)= 
R6(1 54)= 
RG  Cf5T)= 
RG (1 56) * 


0 

0 

0 

a 

o 


RGC158)*  0 
RGtT5D  =  O  ' 
RG (1 60)=  Q 
RGC16D*  0 
RG<162)=  0 
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TABLE  B-l.  MODIFICATIONS  TO  BOUT  AND  RVOUT  (Continued) 


P6(  4) =  PGC  4?»RT0P 
PG(1 0) a  P6(10)*RTOD 
P6<16)a  PS(16)«RT0P 
PG (43) »  PG(43)*RT0D 


PG<52>  a  PG ( 52) *  RTOD 
PG  (  3  3  )  =  PG ( 33) »  RTOD 
PG(53)  *  PG(58)*RT0D 
P6(64)  =  PG(64)«RT0D 


PG (71 ) * 
PG(77)- 


PG (71 ) *RT0D 
~PG~t77) »RT0D 


PG(83) =  PG(83)*RT0D 


PG (94) *  PG(94)*RT0D 

PG(95  )-  PG(957*RTOD - 

PG(96) *  PG(96) *RT0D 


PG(100)=PG(100)*RT0D 

- PGitftl )«PG (1 01 ) oRTOP  — - 

PG(102)aPG(102)*RTOD 

PG<  47  -  PG (  4) - 

RG(10)  =  PG( 1 0) 

gsno  =  P  G  ( 1 6) - 

RG(43)  =  PG(43) 

RG(49)  a  pq(49) - 

RG (5  2)  =  PG(52) 

- PG<55)  -  P3( 55) - - 

RG(53)  =  PS<  58) 

RG(6<»7  a  F^(64) - 


RG (71 )  =  PG ( 71 ) 

- RG (77)  -  PG ( 77) - 

RG (83)  *  pg(83) 

RG ( 83)  a  PG(S8) - 

RG(89)  =  PG( 89) 

PS<90>  =  PGX907 - 

RG(94)  a  PG ( 94 ) 

- RG(95)  a  P9(»5) - 

RG (96)  *  PGC96) 

PG  (97)  a-  PG ( 97) - 

RG(100)a  PG(100) 
PG(t<M>*  PG(10tt 
RG(102)a  PG(102) 
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TABLE 

B- 

-1.  MODIFICATIONS 

PM (54) 

3 

UDACB 

RM(  1) 

s 

PM(  1 ) 

RM(  2) 
~  RM< — 7) 

3 

PM(2) 

'S~ 

r  M(  3) 

RM(  4) 

3 

PM(4) 

-  RM( -57 

3 

PM<53 - 

RM(  6) 

3 

0 

RM<  77 

*  PM  (6) - - 

RM (  8) 

X 

PM(7) 

—  RM-<-  97 

3- 

PM<-8) - 

RM(10) 

3 

PM(9) 

RM(1 1 ) 

= 

PM<tO> - 

RM (1 2 ) 

3 

0 

RM<1 37 

3 

PM<t"1>~ - 

RM (1 4) 

3 

PM(1 2) 

-  RM(15)~ 

PM  (737 - - - 

RM (16) 

= 

PM ( 1 4) 

RM(1 7) 

3 

PM(15) 

RM(1 3) 

3 

0 

RM(1 9) 

X 

PM(1 67 

RM (20) 

3 

0 

RM<21) 

X 

o  - 

RM(22) 

X 

PM(17) 

RM(23) 

3 

PM(t8) 

RM (2  4) 

X 

PM(19) 

RM(25) 

X 

PM( 20) 

RM(26) 

= 

0 

RM(277 

X 

0 

RM(28) 

= 

PM ( 21 ) 

RM(29) 

X 

PM( 22) 

RM(30) 

3 

PM( 23) 

RM(31 ) 

3 

PM( 24) 

RM(32) 

3 

0 

*n  \J  J) 

S 

0 

RM(34) 

3 

PM( 25) 

RM(35> 

3 

PM( 26) 

RM(36) 

3 

PM(27) 

RM(37) 

3 

PM(28) 

RM (38) 

3 

PM( 29) 

RM(39) 

3 

PM ( 30) 

RM (40) 

S 

PM(31) 

RM(41) 

3 

PM(32) 

RM(42) 

3 

0 

RM(43) 

3 

PM(33) 

RM(44) 

X 

PM(34) 

~  RWT457 

~a 

PMC757 - - 

RM (46 ) 

3 

PM ( 36) 

RM(47) 

3 

PM(37) 

RM(48) 

3 

0 

TO  BOUT  AND  RVOUT  (Continued) 


RM(5Q)  »  PM( 39) 
RMCS17  a  PM(40) 
RM(52)  a  PM(41) 
RM(53)  *  PM(42) 
RM(54)  =  0 
RM(557  =  PH(43) 

RM (5  6)  *  PM (44) 

~  ~  R*c?rr  a~pnt4sr~ 
RM(53)  *  f>M<46) 
RMCS97  a  pim<47) 
RN(60)  *  o 
RM7617  a  PWC48) 
RM(62)  a  pm(49) 

- RMT63)  =  o - 

RM(64)  *  o 
RMC65T  *  T>M(507  ' 
RM (66)  a  pm( 51 ) 

_  RMT67T  a  pmr52) 

RM (68)  a  pm( 53) 
RWT697— P n(  747  - 

RM(70)  a  o 
RM(7T)  *  0 
RN(72)  a  o 


APPENDIX  C 


U70  PROGRAM  .CSS  (Command  Substitute  System)  SOURCE  FILES 

The  Perkin-Elmer  3220  computer  on  which  the  U70  program  was  run,  has  a 
development  compiler  which  compiles  rapidly  but  produces  slow-running 
object  code,  and  an  optimizing  compiler  which  compiles  slowly  but  produces 
fast  running  code.  The  development  compiler  was  used  to  compile  the  U70 
program,  since  compile  times  rather  than  execution  times  were  at  a  premium 
during  the  period  of  performance  of  this  task.  As  stated  in  Section  III,B, 
in  order  to  further  reduce  compile  times,  each  subroutine  was  assigned  its 
own  data  file  so  that  one  could  recompile  one  subroutine  at  a  time  rather 
than  having  to  recompile  the  whole  program  each  time  a  change  was  made  in 
one  of  the  subroutines. 

There  are  two  .CSS  commands  that  can  be  used  to  actually  carry  out  a 
compilation:  U70CA  and  U70C.  The  U70CA  command  compiles  all  the  Pershing 
II  source  files,  while  the  U70C  command  is  used  to  compile  an  individual 
subroutine.  The  specific  job  control  sequence  which  initiates  the  U70  com¬ 
pilation  is  as  follows: 

U70  @1,  @2,  @3,  @4,  @5 

@1  Input  source  file 

@2  Object  file  uam'»  (defaults  to  @1.0BJ) 

@3  List  file  name  (defaults  to  @1.LST) 

@4  Start  option(8>  sepa.ated  with  spaces  (defaults  to 
no  options) 

@5  Load  size  for  F7D  (defaults  to  50K  bytes). 

A  cross-reference  listing  can  be  obtained  by  including  the  command  XREF  in 
the  @4  list  of  start  options. 

A  preprocessing  program  called  TRAJ  must  be  run  before  the  U70  program 
is  ready  to  run.  The  TRAJ  program  is  derived  from  equations  resident  in 
the  Pershing  Launch  Computer,  which  is  located  on  board  the  missile.  Given 
the  launch  and  target  latitude,  longitude,  altitude,  and  the  number  of  sta¬ 
ges  (Table  C-l),  it  computes  launch  presets  and  transmits  them  to  the  por¬ 
tion  of  the  program  which  simulates  boost  and  terminal  guidance.  This,  in 
turn,  generates  a  data  file  containing  estimates  of  the  desired  velocities 
and  cutoff  angles,  estimates  of  flight  and  pitch-over  times,  and  target 
offsets  (Table  C-2).  The  U70  program  uses  these  outputs  to  make  its  own 
calculations,  interpolating  the  TRAJ  output  tables. 

Execution  of  the  U70  program  is  initiated  using  the  routines  UST2  or 
UST3.  These  commands  invoke  .CSS  files  containing  other  commands  to  load 
the  task,  make  logical  unit  assignments,  and  start  the  task. 
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The  UST2  routine  is  used  in  order  to  make  a  full  run  (flight  from 
launch  to  impact).  The  Inputs  to  UST2  are  obtained  from  the  TRAJ  programs 
and  are  used  to  build  a  new  data  file  (logical  unit  11).  This  file  is 
used  to  create  a  new  input  tape  in  the  U70  program,  TAPE  11.  The  output 
from  UST2.CSS  is  printed  at  time  Intervals  whose  value  is  input  in 
Constant  271,  or  is  printed  at  any  point  along  the  flight  path  where  a 
discontinuity  appears  (such  as  the  end  of  the  boost  phase,  the  apogee, 
etc.,). 

The  subroutines  in  the  U70  program  are  linked  using  the  .CSS  command 
U70LINK.  This  command  builds  an  overlaid  executable  task.  For  example,  it 
ensures  that  copies  of  all  of  the  subroutines  called  within  a  given 
subroutine  are  present  within  that  subroutine. 

TABLE  0-1.  TARGETING  PROGRAM  INPUTS 


ITRS 

-  210 

Single  Stage 

220 

Two  Stage 

LONT 

m 

Longitude  of  the  Target 

LATDT 

m 

Latitude  of  the  Target 

HT 

m 

Altitude  of  the  Target 

LONL 

m 

Longitude  of  the  Launch  Site 

LATDL 

m 

Latitude  of  the  Launch  Site 

HL 

m 

Altitude  of  the  Launch  Site 

TABLE  02.  GLOSSARY  OF  TARGETING  PROGRAM  OUTPUT  VARIABLES 


U70  TRAJECTORY  OUTPUT _  DEFINITION 


140 

GAMDE  +  DELGAM 

Flight  path  angle  desired  before  pitch-over 

142 

AZD 

Flight  azimuth  measured  from  north 

115 

AZD 

*  Same  value 

114 

LATDL 

Launch  latitude 

116 

LONL 

Launch  longitude 

117 

HL 

Launch  altitude 

203 

DLAT 

Almpoint  offset  latitude 

205 

(Unnamed) 

Aimpolnt  offset  altitude 

206 

LATDT 

Target  latitude 

207 

LONT 

Target  longitude 

208 

HT 

Target  altitude 

219 

GAMDE 

GAMME  desired  at  cut-off 

220 

TCO 

Time  of  booster  cut-off 

302 

TFF 

Elasped  time  of  re-entry  vehicle  freefall 

310 

TGAO 

Elasped  time  before  pitch-over 

312 

VRE 

Required  velocity  for  terminal  phase 

314 

TF 

Boost  guidance  flight  time  estimate 

411 

CHI 

Normalized  velocity  loss 

420 

TAU 

Initial  delta  value  time  impact 

421 

AXPU 

Acceleration  begins  pull-up 

412 

DVP 

Velocity  control  bias 
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U70LINK  is  structured  so  that  it  contains  four  overlays.  The  first 
overlay  includes  a  subroutine  called  MINIT.OBJ,  which  presets  all  variables 
to  zero.  The  second  overlay  includes  ENPUT.OBJ,  which  creates  a  new  input 
data  tape.  The  third  overlay  includes  all  binary  files  appropriate  to  the 
boost  phase  of  simulation.  The  fourth  overlay  includes  all  re-entry- 
vehicle  binary  code. 

In  addition,  U70LINK  creates  a  load  map,  U70-U70LINK.MAP,  that  is  use¬ 
ful  for  debugging  purposes.  Control  indicators  are  input  to  set  flight 
options,  and  program  constants  are  entered  to  control  the  initial  con¬ 
ditions,  multipliers,  limits,  errors,  print  controls,  etc.  Aerodynamic  drag 
and  atmospheric  data  are  input  through  the  tabular  data  portion  of  the 
input  file.  A  list  and  description  of  these  data  inputs  is  given  in 
Reference  1. 

After  a  full  run  has  been  made,  a  restart  can  be  created  using  the 
output  from  UST2.CSS  and  the  UST3.CSS  command.  The  restart  procedure  is 
a  very  useful  tool  in  a  study  where  new  conditions  are  imposed  on  a  missile 
at  a  certain  point  in  flight.  However,  the  restart  could  not  be  used  for 
the  nozzle  deflection  study,  since  a  restart  could  not  be  made  during  the 
boost  phase.  Certain  values  which  appeared  in  THRUST  were  functions  of  time 
and  could  not  be  re-initialized  at  the  arbitrary  time  desired  for  a  restart. 

The  values  which  are  required  to  be  input  into  the  UST3.CSS  data  file  of 
UST3.CSS  are  listed  in  Table  C-3.  Tables  C-4,  C-5,  and  C-6  list  other  U70 
files  used  in  the  performance  of  this  task  for  the  benefit  of  anyone  who 
might  want  to  initiate  U70  runs. 

The  procedure  of  outputting  U70  flight  data  to  tapes  required  that 
modifications  be  made  to  the  subroutines  BOUT  (Boost  Output)  and  RVOUT 
(Re-entry  Vehicle  Output) .  The  desired  output  was  written  to  logical  unit 
14  and  a  .CSS  file  was  created,  with  logical  unit  14  set  to  MAG0,  to  write 
the  data  to  the  tapes. 


55 


rABLE  C-3.  GLOSSARY  OF  RESTART  INPUT  VARIABLES 


V  V  V  41  V  V 

U  <U  U  W  U  u 

cSSace 

*H  *H  *H  *H  *H 

*2  *2  *2  *2  *2  2 

o  o  o  o  o  o 

888888 

^  J2  J2  J2  J3 
O  U  V  U  O  O 

e  0  c  c  e  a 

a  a  a  a  a  a 

4  4  (d  4  Q  4 


c  e  e  a  a  e 


a  a  a  >>  >*  >> 

O  O  O  u  u  u 

•H  *H  *U  14  "H  *H 

u  u  u  y  y  y 

•H  -H  -H  O  O  O 

a  a  4  h  h  h 

o  o  o  u  «  « 

a  5.  a  >  >  > 


U  U  U  U  u  U  4)  0)  <U  -H  60  -H 

acoaoar-i^rxuau 

•  «««««  44t04  OH 

flcflBdcccse^^H<o 

000<J00<4<4B 

a  a.  a.  a  a  a,  si  si  si 

aflBBSBunnyyy 

000000411101000 

OOOUUOr^rH^HSSS 

j2' 


OJ  T3  0>  be 

T3  a  -o  c 
a  u  a  a 


U  W  «  a  41  <0 

4)  «)  a  a  a 
so  »  a  u  a  *a 

ID  ID  D  H  Li  0) 

<LJ  'M  H  bO  «H  M 
O  O  LJ  Q  U  <H 

03  O  H  CD 

u  li  h  h  «a  a) 

3  c  *3 

•H  *H  LJ  4J  4J 

O  O  41  41  4)  4 


u  o 
OB  C« 
oi  a 
fl  a 
oi  a 

a  i 

*4  -rl  4)  H 

x-i  jj  a  i— i 

O  T<  3 

|  UUP. 

u  si 

3  U  01  4 

o  -h  >  a 


<U  <4  60 
<u  >  at 
<u  a>  A 
O  O  <4 
I  0  U  0 
U  <6  rH  O 
3  TJ  4)  -H 
y  u  tJ  u 

•u  §>-<  s 

O  (4  « 


«  4  U  « 

B  o  -h  y 


Ed 

5^52 

H 


< 

Q  Q  W  H  M 

53e3£S£ 


ssi  +  3 

£d 


3 

o 

J 

w 

a  H  £->  H 

+  a  S  i  | 

w  J 


u  3 

Q  O  an- 
Sph4S 
<  H  H  H  < 


TABLE  C-4.  U70  SOURCE  FILES 


ARD59D.FTN 

QUAD. FTN 

BEE. FTN 

RVINIT.FTN 

CONV.FTN 

RVAPLT.FTN 

FASTD.FTN 

RVAERO.FTN 

FASTC.FTN 

SCS.FTN 

SORTIT.FTN 

GAUSS. FTN 

SPEEDY. FTN 

RANDU.FTN 

WRTAPE.FTN 

B EXEC. FTN 

BINIT. FTN 

THRUST. FTN 

BGUID.FTN 

TRIHA.  FTN 

FUNC.FTN 

BOUT. FTN 

VACRAN.FTN 

BAPLT.FTN 

FOLD. FTN 

RATLIM.FTN 

GRAPH. FTN 

BAERO.FTN 

NINIT.FTN 

UST2.CSS 

ENPUT.FTN 

UST3.CSS 

RVEXEC . FTN 

UST4.CSS 

STAT.FTN 

U70LINK.CSS 

RVOUT.FTN 

U70C.CSS 

RVGUID.FTN 
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TABLE  05.  FILES  USED  FOR  DEFLECTION  STUDY 


C0B351.DAT 

- 

30 

seconds. 

.5 

degree 

pitch 

COB352.DAT 

- 

30 

seconds , 

2.0 

degree 

pitch 

C0B353.DAT 

- 

30 

seconds, 

7.6 

degree 

pitch 

C0B354.DAT 

- 

49 

seconds , 

.5 

degree 

pitch 

COB355.DAT 

- 

49 

seconds , 

2.0 

degree 

pitch 

COB356.DAT 

- 

49 

seconds , 

7.6 

degree 

pitch 

COB357.DAT 

- 

30 

seconds , 

.5 

degree 

yaw 

C0B358.DAT 

- 

30 

seconds. 

2.0 

degree 

yaw 

C0B359.DAT 

- 

30 

seconds. 

7.6 

degree 

yaw 

COB360.DAT 

- 

49 

seconds , 

.5 

degree 

yaw 

C0B361.DAT 

- 

49 

seconds. 

2.0 

degree 

yaw 

COB362.DAT 

- 

49 

seconds , 

7.6 

degree 

yaw 

TABLE  06.  FILES  USED  FOR  TRAJECTORY  PROFILE 


C0B370.DAT 

COB371.DAT 

COB372.DAT 

COB373.DAT 

COB374.DAT 

COB375.DAT 

COB376.DAT 

COB377.DAT 

COB378.DAT 


APPENDIX  D 

TRW  PROGRAM  FILE  NAMES  STORED  ON  DISK 
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TABLE  D-l.  SUBROUTINE  FILE  NAMES  STORED  ON  DISK 


TWAERO.FTN 
TWAFC.FTN 
TWALNON.FTN 
TWAMC.FTN 
TWATMENT . FTN 
TWATMEXT . FTN 
TWATMOS.FTN 
TWAUX1 . FTN 
TWAUX2.FTN 
TWBATT.FTN 
TWBURST . FTN 
TWCMDPRC . FTN 
TWCROSS.FTN 
TWDCU.FTN 
TWDCUERR . FTN 
TWDERIV.FTN 
TWDMNDC . FTN 
TWDOT.FTN 
TWDRAND.FTN 
TWEBITS.FTN 
TWFLIGHT.FTN 
TWFSICFU.FTN 
TWFSISA.FTN 
TWFSSCFO . FTN 
TWGRAV . FTN 
TWGEN . FTN 
TWIMSALN . FTN 
TWIMSPLS • FTN 
TSIMSRSV . FTN 
TWIMSTST . FTN 
TWINTERP . FTN 
TWJDM.FTN3 
TWMAGN.FTN 


TWMATRX1 .  FTN 

TWMATRX2.FTN 

TWMAXBT.FTN 

TWNAVON.FTN 

TWNESTG.FTN 

TWPWRON.FTN 

TWRADAR.FTN 

TWRCS.FTN 

TWRCSFM.FTN 

TWRGU.FTN 

TWRKUTTA . FTN 

TWRMMUL.FTN 

TWROTATE . FTN 

TWRVAERO . FTN 

TWRVBAS.FTN 

TWRVPTCH . FTN 

TWRVROLL.FTN 

TWRVYAW.FTN 

TWSAF.FTN 

TWSBITS.FTN 

TWSDPROC . FTN 

TWSSICFU.FTN 

TWSSSCFU.FTN 

TWSTINIT . FTN 

TWTHRUST .  FTN 

TWTRNSF2 . FTN 

TWTRNSF3.FTN 

TWTRNSPO.FTN 

TWURN.FTN 

TWVANES  .FTN 

TWVMP.FTN 

TWWIND . FTN 
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TABLE  D-2.  FILES  CONTAINED  ON  DISC  FOR  USER  PROGRAM  HANDLING 


The  following  files  are  used: 

DAVEF7.CSS  -  Used  to  compile  each  subroutine  file. 

TRWLINK.CSS  -  Used  to  link  the  TRW  6DOF  simulator. 

DCOMMENT.FTN  -  Contains  comments  which  describe  the  function  of  each 
subroutine . 

TRWRV.DAT  -  Contains  re-entry  aerodynamic  data. 

TRWSS.DAT  -  Contains  single  stage  aerodynamic  data. 

TRWTAB.FTN  -  Used  to  read  TRWRV.DAT  and  TRWSS.DAT  and  string  them  into  a 
long  one  dimensional  array.  TRWTAfi.CSS  is  used  to  run  this  program  and 
TRWTAB.DAT  contains  the  output  of  the  TRWTAB.FTN. 

BIG01R.DAT,  BIG01I,  BIG01D  -  Contain  variables  defined  in  the  program  and 
their  storage  locations  in  the  BIG01  arrays,  as  specified  by  EQUIVALENCE 
statements. 


The  following  files  are  used  in  the  $INCLUDE  statements  to  handle 
common  blocks  and  equivalence  statements. 

TRW001.DAT 

TRW002.DAT 

TRW003.DAT 

TRW004.DAT 

TRW005.DAT 

TRW006.DAT 

TRW007.DAT 

TRW008.DAT 

TRW009.DAT 

TRW0010.DAT 

TRW0011.DAT 

TRW0012.DAT 

TRWCON.DAT 

TRWCONST.DAT 
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TABLE  E-l.  SUBROUTINE  NAMES  AND  CALLING  SEQUENCES 


AERO 

acu 

FSISA 

NEWSTG 

RVAERO 

INTERP 

SBITS 

TRNSF2 

AFC 

SBITS 

AMC 

DCUERR 

FSSCFU 

PWRON 

SBITS 

TRNMAT 

AFC 

DCUERR 

TRNSPO 

INTERP 

INTERP 

GRAV 

STINIT 

GRAV 

ALNON 

DERIV 

GRN 

RMMUL 

MATRXI 

VMP 

TRNMAT 

AUX1 

IMSALN 

RADAR 

TRNSPO 

NOZZLE 

SBITS 

TRNSF3 

RMMUL 

THRUST 

RMMUL 

JDM 

TRNMAT 

RCS 

AUX1 

TRNSPO 

AMC 

AERO 

RCSFM 

INTERP 

IMSPLS 

CROSS 

DMNDC 

RMMUL 

ATMENT 

SAF 

RGU 

SDPROC 

IMSRSV 

ATMEXT 

BATT 

RKUTTA 

FSICFU 

IMSTST 

DERIV 

ATMOS 

FSISA 

SBITS 

FSSCFU 

RMMUL 

AUX1 

SSICFU 

INTERP 

MATRXI 

SSSCFU 

ROTATE 

GRAV 

BURST 

JDM 

RMMUL 

SBITS 

RAVERO 

MAGN 

CMDPRC 

MAGN 

RVBAS 

ATMOS 

RVROLL 

DOT 

MATRXI 

RVPTCH 

AUX2 

MAGN 

TRNMAT 

RVYAW 

MAGN 

ROTATE 

MATRX2 

DRAND 

RMMUL 

RVBAS 

RMMUL 

TRNSPO 

INTERP 

EBITS 

BATT 

MATRX2 

RVPTCH 

SBITS 

FLIGHT 

TRNMAT 

INTERP 

RKUTTA 

BITS 

AUX1 

MAXBT 

RVROLL 

AUX2 

INTERP 

BURST 

IMSALN 

NAVON 

IMSPLS 

TRNMAT 

RVYAW 

CMDPRC 

IMSRSV 

TRNSPO 

INTERP 

SBITS 

RGU 

MATRXI 

STINIT 

SAF 

CROSS 

FSICFU 

GRAV 

SBITS 

UNITV 

SBITS 

RMMUL 

SBITS 

SDPROC 

SBITS 

SSICFU 

SBITS 

SSSCFU 

SBITS 

STINIT 

RMMUL 

THRUST 

INTER? 

TRNSF2 

TRNSF3 

TRNSPO 

URN 

VANES 

VMP 

INTERP 

WIND 

INTERP 
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APPENDIX  F 

ORGANIZATION  OF  THE  INTERP( INTERPOLATION)  DATA  TABLES 
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APPENDIX  F 


ORGANIZATION  OF  THE  INTERP< INTERPOLATION)  DATA  TABLES 

An  interpolation  routine  called  INTERP  forms  the  backbone  of  the  TRW 
program's  aerodynamic  simulator,  and  utilizes  a  large  set  of  multidimen¬ 
sional  tables  to  generate  its  results.  The  multidimensional  interpolation 
process  is  straightforward  and  is  described  with  the  aid  of  the  following 
example. 

Suppose  that  we  have  a  table  of  pressures: 


TABLE  F-l.  AERODYNAMIC  PRESSURES 


0  0.1 

0.3  0.5  (bars) 

and  a  table  of  temperatures, 

TABLE  F-2. 

AERODYNAMIC  TEMPERATURES 

1000° 

1300°  1600° 

The  pressure  and  temperature  are  considerd  to  be  independent  variables. 
Associated  with  these  Independent  variables  is  a  table  (Table  F-3)  of  aero¬ 
dynamic  coefficients,  C: 


TABLE  F-3.  AERODYNAMIC  COEFFICIENTS 


0.30028 

0.30057 

0.30081 

0.30117 


0.30021  0.30106 

0.30021 _ 0.30176 

0.30021  0.30268 

0.30022 _ 0.30345 


Now  suppose  that  we  are  given  pressure,  P,  -  0.3306  (bars)  and 
temperature,  T,  -  1116°,  and  we  want  to  interpolate  to  find  the  aerodynamic 
coefficients  that  correspond  to  these  values  of  pressure  and  temperature. 
First,  the  Interpolation  routine  tests  the  input  pressure  value  P,  against 
the  table  of  pressures  and  determines  that  the  input  value  P,  -  0.3306  lies 
between  0.3  and  0.5.  The  location  of  the  upper  number  -  the  "0.5"  -  is 
read  out  and  stored  in  a  variable  which  might  be  labeled  "PINDEX"  for 
Pressure  INDEX.  Since  0.5  is  the  fourth  element  in  this  independent 
variable  table,  the  number  stored  in  PINDEX  is  a  "4",  and  it  is  a  pointer 
to  the  interval  in  the  independent  variable  table  in  which  0.3306  is 
located. 
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Next,  Che  same  operations  are  carried  out  for  Che  temperature,  using 
the  temperature  table.  Since  the  temperature  lies  between  1000s  and  1300°, 
the  location  of  the  1300s  value  -  location  #2  -  Is  stored  In  TINDEX 
(Temperature  INDEX)  to  designate  the  proper  Interpolation  Interval. 

Finally,  the  two  Index  values,  4  and  2,  are  used  to  locate  the  Inter¬ 
val  around  the  4th  row  and  2nd  column  of  the  dependent  variable  table  of 
aerodynamic  coefficients  (Table  F-3).  Then  the  pressure  and  temperature  to 
be  evaluated,  P  -  0.3306  bars  and  T  ■  1116s,  are  used  to  carry  out  the 
actual  interpolation,  in  order  to  obtain  the  interpolated  values  of  the 
aerodynamic  coefficients. 

In  the  listing  of  the  INTERP  subroutine's  data  output,  the  following 
commands  appear: 

BEGIN 

INDEPENDENT_TABLE  OUTPUT  VARIABLE/ 

INPUT  VARIABLE 

DEPENDENT_TABLE  OUTPUT  VARIABLE/ 

INPUT  VARIABLE,  INPUT  VARIABLE 

END. 

At  first,  these  were  thought  to  be  VAX  11/780  commands.  Later,  it 
was  determined  that  the  TRW  program  contains  its  own  data  input  processor 
and  logical  analyzer,  which  reads  and  Interprets  input  data  tapes. 

Including  the  above  statements.  Apparently,  these  table  definitions  are 
used  only  for  printer  output  formatting.  Using  the  example  given  above, 
these  statements  would  take  the  following  form: 

BEGIN 

INDEPENDENTJTABLE  PINDEX/PRESUR 

INDEPENDENTJTABLE  TINDEX/ TEMP 

DEPENDENT_TABLE  AER0C0/PINDEX,  TINDEX 

END 

When  the  results  are  printed  out,  AER0C0,  standing  for  "Aerodynamic 
Coefficients”,  would  be  printed  at  the  top  of  the  column  listing  the  aero¬ 
dynamic  coefficient  values  which  INTERP  generates  by  table  interpolation. 


66 


APPENDIX  G 

FLOWCHART  OF  AN  EARLY  MAIN  PROGRAM 
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MAIM 


Prompt  to 
user  for 
source  file 


Figure  G-l. 


FLowchart  of  Early  MAIN  Program. 
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ONETIM  -  FALSE 


^  SFILE 
+ 

JBLNK(1:40) 


ERRFLG 

* 

10 


ERRFLG  -  0 


IDRNAM  -  SFILE 
SFILE  -  BLNK( 1:40) 
EXTFLG  -  1 

::zz:r~=: 

CALL  REDIDR 


ERRFLG 

* 

0 


d  ;  i  e  ;  i  f 


Figure  G-l.  Flowchart  of  Early  MAIN  Program  (Continued). 


Figure  0-1.  Flowchart  of  Early  MAIN  Program  (Continued). 
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Figure  G-l.  Flowchart  of  Early  MAIN  Program  (Continued). 
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Figure  G-I.  Flowchart  of  Early  MAIN  Program  (Continued). 
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Figure  G-l.  Flowchart  of  Early  MAIN  Program  (Continued). 


APPENDIX  H 

FLOWCHART  OF  THE  SUBROUTINE  SDCTRL 
(SIX-DEGREE-OF-FREEDOM  CONTROL)  PROGRAM 
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SDCTRL 


,SE 


TRUE 


FALSE 


INITFL 


Figure  H-l .  Flowchart  of  SDCTRL  Program 


Figure  H-l •  Flowchart  of  SDCTRL  Program  (Continued) 


Figure  H-l.  Flowchart  of  SDCTRL  Program  (Continued). 
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Figure  H-l.  Flowchart  of  SDCTRL  Program  (Continued) 


m 


Figure  H-l.  Flowchart  of  SDCTRL  Program  (Continued). 
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Figure  H-l.  Flowchart  of  SDCTRL  Program  (Continued). 


APPENDIX  I 

TRW  PROGRAM  OUTPUT  VARIABLES 


TABLE  1-1.  TRW  PROGRAM  OUTPUT  VARIABLES 

AEARTH  5 1  J  0 1  P  (  1  )  EDUATO®  IAL  EART-i  RADIUS 

ALT  31  3C 1 P  C  2 ) ALTTTuC  E 
AST  i IS U 1  R  (  3)ALTITUD?  0  F  TARGET 

DCM  I  3IG01RC  4) INITIAL  CP  OFFSET 

DElTIM  3 1 301 P (  7)INC,REMENT  FPOM  LAST  TAC®  IC  EE  ADDED  SIMTIP 
CEfi  9T50‘1kC  STSTaTE  VmRTAALE  DERIVATIVE'S  ' 

olfsv  uiggirc  2Df:c:t  stage  vane  defl.  angles 
DlSVV  BI30I?r247?V  VANE  OEFL.  ANGLES  ' 

EFTOOf  3I3C1RC  2?) EARTH-FIXED  TO  OESIRZO  TRANSFORMATION  "ATRIX 
EFTCMF  3ISOTRC  57) E AFTH-F IX  ED  TO  MISSILE  T°  An  S  F  OP  ?•’ A  Tl  ON  MATRIX 
E  F  T  0  N  D  sISCIRC  -»C)  EA®TH-FIX  EC  TO  -.OCAL  NED  T®  ANSF  OBLATION  MATRIX 

FRCS?'  3733 1  SCT53  3 C S  TR RUST'YEcTTR- TTj  MI S'S TL7  “FR'AME  — - 

GE  sISCIR (  aG) GRAVITY  VECTOR  IN  EARTH-FIXED  FRAME 

G  m  5  I  301 R  C  62)  ... 

HEAD  31  SOI  R  (  a  6  )  A  Z I  >1 U  T  h  OF  MISSILE  LONGITUDINAL  AXIS  ON  EL 

LAL  alGOTRC  74) L AT  IT  JD  E  OF  LAUNCH 

lAL.A  aI3G1BC  75)  LATITUDE  OF  LAUNCH  IN  ®ADIANS 

LOL  3IG7TRC  7o)L0N’3TtlD'c  "C  F  L A U N CR  ~  . 

LUcR  3 1  3  0 1  R  C  77  )  LON  3ITUDS  C-  LAUNCH  IN  RADIANS 

/DOT  2IJU1SC  73TPAI.N  ENGINE  HASS  FLOW  RATE” 

HFTOEF  3 1  3  0 1  S  (  7':)WIS3Ii.E  TO  EARTH  FIXED  TRANSFORMATION  *ATRIX 
•MJDM  3I3U1RC  2  S )  J  E  T  uA MOMENT  VECTD®  TN  MlSSTLE  FRAME 

iM.RCSr  3I3C1?C  ®1)®CS  MOMENT  VICTOR  IN  -TSSILE  FRAME 
•MTY  — 1I3CTR  T'T41PRC3ULSIV  I  TH  RlTs,p  MOMENT  VECTOR  MIS  SILE  T  RA  ME 
NOZSTA  31 301 S (  97)nOZILE  alVCT  BCI\T  LCCATICN-X  COMPONENT 
P  alSCIPC  STROLL/  =ITCHr  AND  YAW  BCDY  «?ATES-P,  5,  AND  P 

HR  alGOlRCI 32) ATM 3SPHERIC  PRESSURE 

FJI  31  301  R1T  33>  EULER  A  N3LE  3  -  3L  AT  FC  ?M  TO  *1  3  S IL  E“  FR  AM  T 
RE  31 301 fi  <10?)C5lATE  EARTH  RADIUS 

TFRVS  E I  TOT5  Cl  17  >TIM  E  FROM  P  V  3  Z  3  A-R  *  TT  O  N  ' 

I"  iloCIRCI  |R)'A3IC  0  PC  °L  LS  I  V  3  THRUST  VECTOR  IN  MISSILE  F  ®  AME 

TNEw  JlSuT^CI  C2)TTME  AT  END  OF  CURRENT  STEP  ( 3 v 1 TTT  *  D  ELTIV) 

TDL  3ISC1R(123)TIME  OF  LIFT-OFF 

TORVS  UI 3ClRC124)TIMr  0s  RV  SEPARATION 

TJ3I  jIj01SC125)TIME  OF  SECOND  STAGE  IGNITION 

U  -UIGC1»T13t) VEL OCITY • IK  M I  STILE  FRA^E  . 

3IGC1 R C129)c4PT*  ROTATION  3ATE 
XCM  STA  313017(1 50)  MISSILE'  TRANSIENT  CM  LOCATtCN 
XE  alGul R (1 31 ) °0S ITICN  I'  EACTH-piXED  FRA«E 

XLE  r  I  SCI  R  Cl  34  )  L  AL  NC’H  3OSITION  VECTOR  IN  EARTH-FIXED  FRAME 

aTE  j.1  3C  1  ?  (1  37)  TAR  SET  '’OSITIDN  VECTO3  IN  EARTH-FIXED  ®RA*E 

AuFDOT  SIOCrPTIAO  )T-Mr-DE’IVaTTVE  0--AL7HAM- -  '  “ 

A  l  H  pi  A v  2I3j1cC1“1)0ITCH  PLANE  COMPONENT  0®  ANGLE  Or  ATT  AC< 

ALPHAT  TTGOTF (142) TOTAL  ANGLE  CF  ATTACK 

A1CDLX  a  I  3  0 1 R (1 43)AISD.  MOMENT  COEFFICIENT  TP  ANSF  E°  DISTANCE 
A1CSTA  3I3C17C1-.A)  AERO.  70MENT  CCEFFICI  ENT  REFERENCE  LOCATION 
AJL  3ISu1F(1i5) ALTITUDE  OF  LAUNCH 

*0L  ■-  31  SCI  RC145  )  ALTITUDE  0*  LAUNCH  ‘  " —  "  “ 

b?.-r  Jl  SCI  p  Cl  46)  AE®0.  3  EF  ;?  £NCE  A  ’E  A 

-ETAV  3ISL1PC1V7)YAW  PLANE  COMPONENT  OF  SN3LE  OF  ATTACK 
_  a  T  D  C  T  313  71 P (1 42 ) TIM E  DERIVATIVE  OF  3FT  *  * 
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TABLE  1-1.  TRW  PROGRAM  OUTPUT  VARIABLES  (Continued) 

CAT  31 SCI  R  Cl 4? ) TOT  AL  4330.  AXIAL  FORCE  COEFFICIENT 

Cf*  315010(1  ;C)cA5IC  N0P1AL  FORCE  COEFFICIENT 

CN - mClFCt  TC7D~ATAi*ET<ER  FO1*  ATCTLTFATT  ON- LINT! - 

CNSTAT  31  SCI  S  Cl  51  )  TjT  4L  STATIC  NORMAL  FORCE  COEFFICIENT 
CNT  TIGC1>^CT5? JTOT.AL  AERO.  NORMAL  “CTC“~  COEFFICIENT 
LPMT  9I301R C152JT0TAL  AERO.  PITCHING  MOMENT  COEFFICIENT 
C  RMT  31 GTOTR  CT5  4  TTOTRL  ~  ATFC.  RUtrLTNC  “MTTMTNT  CUFF  FITTETIT 

CS  alGOI  R  (1  55  )  SPE  EO  OF  SOUND 

CY  "  3TGCTR C1;6  OEASTC^STtrr’nJRCF  OUr'TTTTEKT'  -  _ 

CYNT  31  GC1  R  C157)TOT  AL  AERO.  YAWING  FOMENT  COEFFICIENT 
C Y STA T*  SI  301~R  tT5 3 0 T  JT AT.  STATIC  STOE  FORCE  CTEFTICTENT 
CYT  3I3C1R(159)T0TAL  AERO.  SIDE  FORCE  COEFFICIENT 
OErRV~  31  SOI  F  Cl  St  DRV  VANE  OECrECTTCN  OS  IV~  T  N  TNTE  RTOL  AT  ION 
DELTAP  3IGClRC1o1)E3LI VALENT  ROLL  CONTROL  VANE  DEFL. 

DELTAS  sIGUlRCIoDtiLIVAuENT  HTIC-I  CONTROL  VANE  DEFL. 

DELTA*  31  SCI R (1 32 ) 50L IVALENT  YA-  CCNT°CL  VANE  DEFL. 

LLFSVP  sIG jIRCI 54) CIR3T  STAGE  ROLL  CONTROL  VANE  D=FL. 

DXDRL  iIjC1?(16S)TOTAL  AERO.  MOMENT  CCEC.  TFANSFER  DISTANCE 
DYNPR  a  I SCI R (1 6c)DYNAMlC  PRESSURE 
i-I  ?T  E  31 301 R  (1  37) ELL  IPICITY  OF  E*RTH 

E  1  S  A  is  =  I  GoT  R  (  T  5  3)  3  I A  S  ERRORS  IN  TOTAL  A  £  R  0  .  C  C  E  F  F I  C  I  c  N  T  S 
C.PSAF  3I301AC174) FRACTIONAL  ERRORS  IN  TOTAL  AERO.  COEFFICIENTS 
EPSCY3  3I3C1R  (1  2  j  )  £  I A  S  EF.AG°  IN  LDNOITLDI  N  A  L  CM  POSITION 
E  PS  IS  3I3C1  ?  Cl )  EIA  3  ERRORS  IN  MOMENTS  »ND  PRODUCTS  OF  INERTIA 
-  =  3ir  31  SOI  PCI  3  7  J  ERRC’O  IN  MOMENTS  AND  PRODUCTS  C'F  INERTIA 
cPSjls  =  I3C1  ?  Cl  R?)W'AS  S  EIAS  E’RSS 

EPS  V~E  3  ~3T  3C1  ATI  ?4D  EfI  A  j  ERROR  IN  v0MEM  CCEF.  OLE  TO  VANE  “DEFL. 
tPSVEr  BIjCIRCI^TJE^oOR  IN  mc*cST  COEF.  DUE  TO  VA.NE  DEFl. 

F  AM  s  I  SGI  F  C  EOC  )  A  £R  0  .  FOPCE  V’CT"CR  IN  MISSILE  FRA 'IE 

f*M  s I  3  0 1 P CIDiJMIcSIwE  WEIGHT  VECTOR  IN  1I3SILE  rPAM£ 

ITEn  "31  SGI  A  CZZf)  INERTIA  TEn.'R 

ITENCT  ilSul A C215)TIME  RATE  CF  CHANGE  C?_mCMFNTS  OF  INERTIA 
LAE  r  JISUIRTOI5) AE’0.  REFERENCE  LENGTH 
MALH  31  SCI  A  CZ  20  MAC  H  NUMBER 

MAM  31301 F CZ21 ) A EFO.  MOMENT-  A  3  O'JT  MISSILE  AXES 
MASS  EI3C1 R CZZ4J =0C ST  *HSS  ESTIMATE 
•’ASS  a  I  3)1  R  C  2  Z  4  )  V  EH  I C  LE  WA3S 

i'  A  S  S  G  3  3 1  3C1  3  C  Z  2  5  )  M  A  0  S  EX0ELlEC  SY  SAG  GENERATOR _ 

OG'c  Jl'SCI  R  C2  !i)  GAS ‘G  E'*!  =  R AT  OR  MhSS  -LOW  RaTE 
r  H I A  P  EIoUlR  CZ2o)A£FJ.  pJLL  ANGLE  MEASURED  IN  A/P  FRAME 
PRIOR  aluClR C22=) AERO.  ROLL  ANGLE 

P.-tIHV  SI  30 1RC23D  MODIFIED  a£RC.  ROLL  AN&.-E  FOR  INTERPOLATION 
r  HO  EISCI  R  (231  )  ATMOS=»h£R  IC  DENSITY 
RuOZV  _  3IGQ1F  CZJDLREF  /  C2.C  *  _V  AT) 

*  '*  “  3  I  3 OTR  C  c  5  2  )  A  c  Y  N  G  L D  S  N L  *  r.  r. « 

HNTHNS  31  SCI R  C  2  3  4 ) r  CL  ND  ARY  LAYER  TRANSITION  PEYNCLDS  NUM3ER 
SUMFM  i  TSC  TP  (  2  2  5  )  S  Jv  M  A  T I  ON  CF'  FORCES  IN  FISSILE  r  FA*  E 
i  jM-MM  G 2  SLIP  (25;  )  sUvwATION  C1  IN  MISSILE  FRAME 

T  eL  ?Ij01a(2C1)TIvE  FRO*  LIFT-OFF 

t  :  s  i  =:  sui  -  c:  '*i )  t  rv  £  =  =  q,v  :£cond  stage  ignition 
TOLD  ,rs^1  A  (2  •»  7 ) 
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TABLE  1-1.  TRW  PROGRAM  OUTPUT  VARIABLES  (Continued) 


-£F  iIoLlD(.44)'vI?3IL£  INERTIAL  VEL.  *.9.  T.  E  A  9T  FIXED  FRAPS 
VAT  713G19  CZ'47)MA;TNr’UDE  OF  M^Slt:  VELOCITY  W.P.T.  AIR  " 

V  s  2I3C1R(243)PISSILS  E A R TH- R EL AT  I V E  VELOCITY 

VE - TI  GOT  9 C243T EARTH"* £L ATTV E  VELOTTTT  VECTOR  - 

VI  3 I3C1R(249) MISSILE  INERTIAL  VELOCITY 

VT3r - 3 1 131  ST  2  jCTHETTin  FTT"C'F"  VTTCTSTTY  . 

VuAT  _iIGu19(251)VEHICLE  LATITUDE 

VLOfJ 21 371 9 (23  Z )  V  EPIC  lE  U3V3TTUD7  '"  “  " 

M I  is*  t>  d  21  301  ?  (253)WIND  VECTOR  IN  EARTH- FIXED  F9  APS 
XEVE  "sIGOTrt  CcTI  >  c  ATTh  9  ELAT  IVT  TElCC  ITT  COMPONENTS 

ASXIT  iI301fi(229)N02ILS  EXIT  AREA 

DTP"" - ITTDTR  CTT CT F  IUTf  TTTG e  "ATJZTCTT^TT CT"C  C’MM  WTT 

5 1  Y  =12019(211 ) FIRST  STAS-  NOZZLE  Y  A  V  COMPAND 

DZP "313719  C7T7V5ECDND  STAGE"  NOZZLE  PITCAT'CTMA  NO . 

CLY  2I3C19(313)SEC0ND  STAGE  NOZZLE  YAW  COPMAND 

53T  =  17019  (314TTV  YANc  ~T"C  OPMAND  ~  ""  ""  *  " - 

D32  _  3IS01P(215)RV  VANS  2_C0P"1AND__ 

07  3'  7I7oT9  C3Tt )  ?  V"V A  \ E~3  fTK'^D  ''  

554  21  Sul  9 (317 ) 9V  VANS  4  COPLAND 

DC  j  “ 3 1 GGIfi  C7T3TCTS  =  C  TION'C  OS  I.N’ES  CF  TnE"3CS  'yHTU7T  VECTORS 

5i_FSVD  E I  301 9 ( 342 ) F 19  ST  STA5£  VA NE  INSULA3  RATES 

OLNP  ~  3I301R  04 iTNOlILE  ^TTCh  AV3LE 

DlINPO  3I3G1  9  (347  )N0:ZLE  °ITC-t  A  N  3  ’J  L  A  °  RATE 

DLN'Y - .VT  3  GTE  CT  47771777  L' = ~  Y  AO.  TN7LZ  - 

DLNYD  21  =  01  9  C349  )N  02  ZLE  Y  A  *  ANGULAR  RATE 

DL<WvD  31  GOi  P  (750)  R  V  If-  NE  "A  NC-'JL  AR  RATES  . ' 

D  R  C  IIGul 3  (' 54D FIRST  STA3E  VANE  3CLL  COPLAND 

Df-CS  ■  ilSOTS  C3550RCT7  ET  TIME  DELAY"'  -  - 

EPST.jS  31  3C1  0  (  3  3  C  >THPUSt  3IAS  ERROR 

=  3  S  T  ,4  C  3I5C1R  (357)  FRACTIONAL  ERRORS  IN  R  C  3  TH-»UST  MAGNITUDE 
.  3  S  T  3  F  21301  9  (  ’  =  5  )  FRA  CTIONAL  ERROR  IN  TH9JST 

MC3  '  'I3C1R Cj5c)\CPINAL  MAGNITUDE  CF  RCS  THRUST  FOR  ITH  JET 

c  ■'’0  -I  SCI  F  (374;  GRAVITY  PARAMETER 

J  2  -/IjC19(377)S£C3ND  ORDER  ZONAL  -AiRMONIC 

j  3  3I3C15(3"’e)THI3D  C3DE3  ZON  AL_  H  »  9  HO  N  I  C 

j*  =  1  SCI  9 (577) FOURTH  OFDcR  23  n'AL~  H  A  RP  ON  I C 

LAPinZ  3iov.1F(T7-) 

i.  j  D  =13019(379)04  A  RACTEPISTIC  JET  0  A  M  °  t  N  i  lENGTH  POSITION 
XAj  j  I ' 3  C 1  9  (  2  2/M  SPENT  ap*  V5C  TD3S  OF  T  H=.  RCS  JETS  (P.  JETS) 

N  L  I  M  F  3  ilC-ClP  (4.3c)  FIP  ST  STARE  NOZZLE  DErL.  ANoLE  LIMIT 
NuIPiOS  513219(407)  SECOND  STAGE  NOZZLE  0  E  F  L  .  A  N6  Le  j-  IP  IT 
riil  2  2l3JlRC4ji) 

r  3L  V  =I3C1R(4J9)SEA  LEVEL  F3E$SJ»E 

•CSOLX  =1301  3  (41C)9CS  v  0 v  £  N  T  COMPONENT  AL  ON  :  X  AXIS 

*  £  3  ST  A  313-19(411  )  3C«  LOCATION-*  COMPONENT 

*nc:  jIjOI  r  («.2S)t:w  e  at  w-:ch  l-.st  nozzle  angles  computed 

T  V  A  ,» £  S  Si  3J1  9  (4  21  )  TIPS  AT  W^TC*  lASt  VANE  ANGLES  COMPUTED 
Vt-IMf  S'  si  571  4  (422)  fls  St  STAGE  van  e  defl.  an  sl  =  limit 
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TABLE  1-1.  TRW  PROGRAM  OUTPUT  VARIABLES  (Continued) 


Vi-IrfRV  jI3U1R(42Z)'V  VANE  DGFL.  ANGi-E  LIMIT 

*  F  3  V  jIGuIFULOFIRST  STAGE  V  AN  £  NATURAL  “RECJEwCY 
*NP  3I3ClP(4c5  )N0ZZL£  =ITC^  NATURAL  FREQUENCY 

•  NY  '  31 3C1 R C426) NOZZLE  YAW  NATURAL  FREQJENCY 

w  R  V  V  JI3C1  R  (427)  RV  VANS  NATURAL  FR£0UcNCY 

XN07  3T57JTPT42in  IN"!?! AL  NOZZLE  POSITION  WITH  RESECT  TO  CM 
ZETANP  3I3G1 fi (471 )NOZ2L£  °ITCH  DAMPING  FACTOR 
ZETANY  32GOlRU32)NOZZLS  YAW  DAMPING  FACTOR" 

_  ^  T  F  3  V  alGCI P(43?) FIRST  STAG:  VANE  DAMPING  FACTOR 
:  z  TKW~S  IS  31  R  C5 34'  >  ?  V  VANE  DAMPING  FACT  OR 
A ; I  A  S  3I301R ( 4  5  3  )  A  C  C  EL  IP  CM  ET  ER  STASES 
A  3P  3 ToTT?  rooOORAVimiONAL  ACCELERATION  VECTOR 


/*  jP0LJ_  =  1  301  R  (45'A  )  PREVIOUS  GRAVITATIONAL  ACCELERATION  VECTOR 
'  '  ilGLI  5  (4577*11  «UT"H  TYA  WTRESClV  ER  OFFSET 
AGO  sIGC'l  P  (4©2)  AZIyLTH  (YAWORcSCLVER  OFFSET 
DDYT  ~  SHGTTCToTJY  TT  X  NCVORfHCGOVALITY  P=E3TT 
•OOYX  S 1  3  C- 1  R  C  -*  t>  5  0  Y  TO  X  NONC  “ThCSCNALITY  “RESET 

DELYX  Tl  GO  T  S  C 4o 5T7.TN CRTifO S  0 N  AUTTY  i ETWTTN~  Y~" A ND~  X’  AO  C  ETTRCMETE  R S 
OiLIX  il 301 3 U;7 ) N ON  ORTHOGONALITY  3  E  T  W  £  E  N  Z  AND  X  ACCELEROMETERS 
DELZY"  TimTC45E)MONOFTHOGOKALITY  3  E  T  *  S  E  N  Z  AND  Y  ACCELEROMETERS 
OF  3lGG1?(43<5)3YRA  RESTRAINT  DPIFTS-X,  Y/  AND  Z 

Cl  ’  3IG0TFT47Z)  ZNF'JT  yAS3  -JN3ALANrCE  DRIFTS -X,  Y /  AND  Z 
DOZ  sISGI.l  (4?5)Z  GYPC  ACCELERATION  ALONG  X  CLUSTER 
DO  lIGaTRTC7trrjT.RDT  AAS3  TJ’IE'ALANOE  D^IFTS-X/  y;~  ANO  Z 
EIM3  313 OIF (47*)cwL?R  A NG l E 3- TRUE  IMS  P.ATFCR*  TO  PLATFORM 

:  =  sgs  sizar*  C4  32  rn  actional  esrc9  accelerometer  bias  estimate 

:r5Cij  3  I3C1  R  (4  -  5  )  Gif' 3AL  ANGLE  RINSES  FO9  COSINE  RESCLVEPS 
E 3 S  CFG  3I3C1  F  (*.;?0  FRACTIONAL  "ER  R05  IN  IMS  COSINE  RrSCLVER  ME*  3L9E 
. r  3  r  D  -!l  301  R  U?1  )  ERROR  IN  PLATFORM  FIXED  DRIFT  ESTIYiTES 
E  =  S  30  STGIT?  C47  4TF"C~ATTTCN’A'C~  ERROR"  TN  1^3  ERTIT I  V E"  DRIFT  “ARTELS 
: J  3  N  S  F  31 3C1 R <4R?) ERROR  IN  NcG.  A CC EL ER DM ET E»  SCALE  FACTOR 
E  =  S  3  R  T  31  JC1R(50C)E  =  P05  IN  AC  C  cLE'RCMITc  R  NCN3RTH03CNALTIY 
-  3  S  P3  F  ilGCI 9 (503) ERFOR  IN  =  C3.  A  C  C  EL  E  ROM  ET  =  P  SCAlE  FACTOR 
:  =  3333  31 3C1?  CrOo)  GIMAL  ANGLE  EIA3E3  FOP  SINE  RESOLVERS 
:PiiFG  Jl 30 1 R (5  09  ) FRACTIONAL  E  R  5  C  R  IN  IMS  SINE  RESOLVER  MEA$L=E 
E  =  STD  ErGJT3TS'TZTERR'3“"  Tn  TCRCjER  $U£w‘‘F"5TEr  - 


1  a  s  T I M  -I j01R(516>I1S  TEST  TIMER 

ItTOIF  RI3C1RCS17)  INERTIAL  REcERENCE  TC  I  N  E  RT  I A Ea  ATM  TP  AN  S  . 
IiTOPF  jIG'.T  s  (SZs)  IN=  RTIAL  P£Ft<c\CE  TC  PLATFORM  TR  AN  S .  MATRIX 
\l  jISGI  9(535)  ACCEL^OMETER  LOW  GAIN  SCALE  ERRCRS-XYZ 

a.  I H  i  I  3  G 1  5  CjJi)  ACCELEROMET  ER  HI3M  GAIN  SCALE  =9P0RS-XYZ 


<*L3F  sIGLl  *  (5**2)  NOMINAL  NEGA7IV1  LOW  GAIN  A  C  C  EL  E?  OM  ET  E  R  SCALE 
<0  3I3C1RC54:0»CCELEROM;TEO  LOW  CAIN  31  AS ES-X,Y, AND  I 

<  v  ri  =:  3C1  R  (;4c  )  ACC  ELEROMET  =  fi  hIgd  GAIN  PIASES-X,  Y  r  AND  Z 
<=^SF—  31  JC  Tc  (T4  -  J  NO^TFrAL  P  OrTTVT  G3TTFT  GAT  N"~  TCI  EL  E  9  C  v  E  T  E  R  SC  5L” 
K.PL3F  jIjj13(35C>  NOMINAL  “OSITIVE  LC»  GAIN  ACC.  SCALE  FACTOR 

xi  j  I  3  “•  1 *(351 )aCC£lE°OvET  £  s  L;«  GAIN  ZYyM=T=Y  f5RD»>-YYI 

.1  GC  1  =  ( :  34)  ACC  EL  E  =  OvET  E®  HI3P  GAIN  !Ym;Tpy  =RRO=S-XYZ  • 

<  *  ilGjIPCEEZJTYR)  T  C  R  C u E =  SCALE  :a;top  =  A ’  10 

W-.T  ?  I  G  J 1  r  C  =  3.)  LATITUDE  OF  tuvt 
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TABLE  1-1.  TRW  PROGRAM  OUTPUT  VARIABLES  (Continued) 

LOT  5I30lR<5a1 )L0N3ITUD£  Cf  T  A  SO  ET 

M  F  T  0 1 R  iI3Q1R(5  i  2JJ*  ISSILE  TO  INERT  hl  REFERENCE  TRANS .  MATRIX 
MFTOND  3ISL13C571)!»iS3IL£  TO  LOCAL  NED  TRANS.  *  AT  3 1 X 
MFTOPF  5 1  3  C 1  P  (530)  FISSILE  TO  PLATFORM  TRANS.  MATRIX 
PfTOIR  3I3C1 R (5 >6) PLATFORM  tc  INERTIAL  RE  FERENC c  TRAN  S .  MATRIX 
PRO  dISO’fi Ccj5)PI7CH  RESOLVER  OFFSET 

PSIDTP  2ISCTR(63a)PREVT0US  VALUE  OF  PRE-LAUNCH  EULER  ANGLE  RATES 
RCPCH  3 1  3  0 1  F (o39) COS INE  IMS  “ITCH 
ACRJl  3 1 G  C  T  F  (  61  G  )  C  0  S  I  f«  E  IMS*  3CLL  " 

RCYAt.  uIGOIRCelDCOSINE  IMS  YA  W 

ALL  31601 R  Cal  2) LAUNCH  LOCAL  EARTH  RADIUS  . . 

RET  _  3ISQ1R(613>TAPGST  LOCAL  EARTH  RADIUS 

RnG  s 1301 “ 161 4) HIG  H  SAIN  ACC  ELE SOMETH RD SI FT*X  AND  Y 

? lG _ £1301 S (6 It  PLOW  G  A IN  ACCGLtRCHET ER  OFIFT-X  AND  Y _ 

RAG  £  I  30 1 R ( 1 1 o ) R 3L L  RESOLVER  OFFSET 

RiPCH  5lGClR<a19)SINS  IMS  “ITCH 
HSROL  31301  R  (tiG)SINE  IMS'  PClL  ' 

RSYAt,  51301  fi  (o  21  )  S  IN  E  IMS  YAW 

SI3DLV  3l30lRCt2zTSTANDARD  DEVIATION  OF  DiLTA  V  SENSED  #Y  1*$ 

31 3 RD  3IG01R C  e  23  >  D  EV  I  ATI  N<?  PLATFORM  DRIFT  CRDM  RANDOM  ERRORS 

si'3R  uV  BlTiTTR'CaTf)  Sf  AND  ARD  DEVIATION  OF  IMS  RESOLVER  VE  A  S  JR  E 

Vi R  a  I  301  R  (626)yISSILE  VELOCITY  IN  INERTIAL  PEFgP^NCE  FRAME 

V P  "  5I3v.1fi(c£9),viTS3lLE  VELOCITY  IN  PLATFORM  FRAME 

VPCLD  EI3C1 R (t  32JPRE VIOLS  VELOCITY  VECTOR  SENSED  =Y  IMS 

AlTP  £T3UT'R'<fc'35>  PREVIOUS  ALT  FROM  C  CV  ERROR  MOD  El 

A  r*h  A  j  ZJ01  B  (  6  2  S  J  T_R  fl  \  Sf  $  R  FUNCTION  -ROLL  ANTENNA  ANGLE 

APS  A  YI  3‘T1  PToaTTTR'A^sFER"  FUNCTION  -VTw  ANTENNA  ANGlT 

^THA  iI3C1R(t3E)T’ANS.  FUNCTION  -’ITCH  ANTENNA  ANGLE 

aPE  '  si  jOTR  CV3$  >  EAST  ERROR  »A>AMET£R  F"R3m  0"cO  ERROR  MODEL 

EPHA  3I301P(64C)TRANS.  FUNCTION  -»CLL  ANTENNA  ANGLE 

3PN - JIGUTSCeATONTRTH  f SR 0 R'FAR  AM  ;T1TR~'F7(C'M'  DCJ  ~ETYRCT  'YODEL 

rPSA  iI331R(tA2jTRAN3.  FUNCTION  -YA W  ANTENNA  ANGLE 

jTKT - ST  JUT  RTo  4TTTR  AXT7~  FU  mC  T 1 0  N— -"PTTC" H  TY\¥S  AN3L  t - 

CR  623ii1«(644)  UPDATED  INcRTlAL  ?CSItION-N,  E,  AND  D(PI) 

D - ArJUTRTo 47)  ERR  OR  'IN  '  3CTITTO"N“N ORTH/  E  A  JT*  4MD~  'DOWN  “  " 

DTACL  3I3C1R(6SC)TIME  SINCE  START  OF  DCU  ACSUISITION  ATTEMPT 
2PS75R'  3T3TTfr  CaaTTMTSTIuE  'JOD  Y  TATT  'MEASURcMT’IT  “FITS  EPTCR 
EPSCER  Jl30lR(c54)Pj  GI*5AL  EIA3E3  FC  R  COSINE  RESOLVERS 
a'PTC'rT'aTJCT rCZZ71  F  T  A  CT ITNAl "  T^TCTrTW  FCT'CCUI NETE'SCl  VTR~  *F*SJ»5 
LPSrGR  3I 3C 1 R (tsC ) F RA CTIONAL  ERRC’  IN  MISSILE  a  OD  Y  RATE  MEASLRE 
cTJJGR  31 331  S  Ct  53TP  j  01  A*. ~  A  NGLe' is 1.*  STS~  FC  R  “SI  N£  RESOLVERS 

:  ?  $  3  F  R  £10 Cl P (ait) eR ACTIONAL  ERROR  IN  RU  SINE  RESOLVER  MEASURE 

c ? E  '  31 301  ?  ( 169  )  DCL  Ts’ROR  c5ST  ‘  ~  '  . . ‘ 

t.RN  3:3ClR(a7G)CCL  ERROR  N0°'u 

n’TiA  71 JC  1  RTT^A  >  TR  A  NrST  E  R“  ?l NCTTTN  ’C“UI« STAN F- RCLL  ANTI  ADGIE  “ 

N  ?  3  A  3I3C1R(t75)TRANSF£R  FUNCTION  CONOTANT-Y#W  ANT.  ANGLE 

<T  H  4  jl  3C 1 R  C  5Tt1TR  A  NS  F  t  R  FUNCTION  C  0  N  ST  AN  T- *  I T  C’H  ANT.’  ANGLE 

v;  £I3LlR(a7i)M4TCM  *3Um'_ItY 

M  jlGGISCiiDFTSJT  STAGE  ROL.  RATE 

r  i  £2  jl31  R  (t  al  )  SEC  ONO  STAGE  3cuL  ®AT£ 

PS  ~T1  JT1R(J32)£IM-ls  STAGE  ’CUl  RATE  ~ 
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TABLE  1-1.  TRW  PROGRAM  OUTPUT  VARIABLES  (Continued) 

P4  6I3G1R(633)RV  EP/A3/SE  POLL  °ATE 

? 5  OIuGIRCs 34)RV~EP  RCLu  RATS 

Pa  31301  P  (oi5)  °V  A5/S3  »CLL  RATI 

FriC  0l30TR(o35)SAU  ROLL  0«TVT~  “  ' '"  ""  '" 

PSC  31 3d  S  (a  57  )  3  AU  YAW  03IVE 

Si - 313rnrrB3eTFIPTT"STAffE"PTTC‘H-RAT5‘ -  - 

SO  3  1 30 1  S  (  o  3  9  )  SECOND  STAC- 1  PITCH  PATS 

S3  3I30rPC‘9C)FTK3LE  STAGE-  8TTCh"RATS 

U4  31  331  R  (  e  91  )  RV  SP/A9/SS  PITCH  R»TS 

3  5  '•  3 1 3  C 1  R  C  5^2  7  TF  FITCH  RATS'-  ~ . 

So  51  301  R  (693  )  3V  Ac  /S  3  PITCH  RATE 

.11  3I301R(694)FIRST  STA3E  YAW  RATS 

A 2  3I3GlR(a95)SEC0ND  STAGS  YAW  RATE 

R3  “3I3C1  R  Ct76)  3IN3LE  STAGS  YAW  RATS  '  ~--  - 

R*  3I301R(697)SV  IP/A3/5S  YAW  PATS 

R 3  3I3GTRC69S)RV  IP  YAW  RATS- 

Ho  3I3C1R (c 99) RV  A3/S3  Yiv  RATE 

PA  ■oT3'Jm70CTULSR  A  vGL:  S-'ANTaNATA'  TC  MaITSILE  FEATfE" 

-iCAPH  iIG01R(70:>C02IM  .ANTENNA  PITCH 
RCARL  3I3C1?(714)C3$lNE  A\t^NI  ROLL  " 

R  CA  Yu  jI3ClR(7Ji)C3SIM  ANTeNM  YA* 

R3APK  '.I3C1RC7:)?5$INS  ANTENNA  PITCH  "" 

R3ARL  313  G1  R  (70j)SIM  ANTENNA  RCLL 

RjAYw  '  3TGC1  <r(7CV)  SITE  'ANTS'IW  YA*  ‘  ‘  ' 

oIGRGR  3I301R  (71 C ) STAND ARD  DEVIATION  OF  "ISCILE  20DY  RATE  MEASURE 
SloP-A  3 1 301, E  (711) STANDARD  DEVIATION  OF  RJ  RESOLVSR  MEASUREMENTS 
THC  3I30lF(7i:>$AL  PITCH  CAIVE 

T-.XAC3  3I301.9(713)HAXIF,UW  TIPS  PC?  AC  3LIS ITTCN- 3T“T)  CL  TN  ALTITUDE 
T1U0N  31301?  (7l4)TIP£  AT  WrilC.H  RADAR  UNIT  TURNED  CN 
CAPON  '3 1  3  G1 UT7TT)  6  A  S'!  C-  A<m~ Ff*ZT  C  05FF:'C1? YT*'C  POVE  V  ON ) 

C=T  3I3G1R(716)SASIC  PITCHIN  3  FOMENT  COEFFICIENT 

C  ?y.i  51301 P  (717)  PIT  CH  IN  3  mows  NT  3  DAMPING  DERIVATIVE" 

CRN  9I30lR(71i)3ASIC  “CULING  POTENT  CCEFciCIEN'T 

CRM®’  i;30TR  (719)5ASIC  ROLLING  POP  ENT  DA^PTN'S  DERIVATIVE 
CYH  3 1  SC  1 R (720) 9 AS  TC  YAWIN3  MOMENT  COEFFICIENT 
CYMR  — 51 37J1F(72T)Y  AW  I7J3  “MOMENT  ?  D  A  F  ®  I N  G'TTR  rV"®TIVE~ 

CCAAE  3I3C1R(722)INCREM£NTAI  CA  DUE  TO  ASP0ELA3.  EFFECTS 
DCADP  3TGC1  FT7  23  )  INC  RE "!£ NT *L "  C  A  DUE  TC  DELTA?  VANS  CEFL. 
uCADfi  )I3G1R(724)INCRSMENTAL  CA  DUE  TO  DELTAS  V  4  Nc  D  E  Fi_. 

OCADfi  a *, 3 C 1  R  (7 Z3 )  I N C  R S7S NT'A L  CA  ’DUE  TC  DELTA-?  VANE  D?*L. 

OCAPCF  Ol301«(72e) ADDITIVE  J  A  SE  PRESCURE  AXIAL  FO»CE  COEF. 

DC  AS  F  5T30TRC727>3<TN  FRICTION  *OCAL  FO.RCT  COIFFITIENT *  - 
DCAVF  jI30  1c(723)lNCsS'1cNTAL  C»  DUE  TC  VANE  MISALIGNMENT 
DCN*E  31301 R (775) INCREMENT  AL  CN  DUE  TC  A  E 5  0 EL A".  EFFECTS 
*  C  *«  D  P  -I  301, R  (7  5  C  )  I  NCR?  MENTAL  CN  DUE  TC  DSLTAP  VANE  DIFL. 

DCNDC  5I3C1?(731)  INCREMENTAL  CN  SUE  TC  DELTAS  T/ANE  TEFL'S 
OCNDW  31301. 3  (732)1  ^ C5EW£NTAL  C-N  DUE  TC  DELTAR  VANE  DEFL. 

0  CNVV  TI3C1® C73I) INC RFMENTAL  CN  DUE  TC  VANE  MI  3  A  LI 3KME NT 
DCP'IA:  31301  R  (734)  I  NC ’ETcNTlL  CPF  DUS  TD  AE°CEL*3.  EFFECTS 
OCPMDP  51301” (7T3)INCRE  NTAL  C3*  OJt  TO  DELTAP  VANE  DEFL. 

DCPMCJ  51301 e  (7-It)  INCR£MFNT  A(.  C5?  DJc  T  0  DElTAG  VANE  DEFL. 
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TABLE  1-1.  TRW  PROGRAM  OUTPUT  VARIABLES  (Continued) 

>CPMLR  3IG01fi(7I7DINCRE'<SNTAL  C?M  SUE  TC  DELTA®  VAN;  DEFL. 

DCPVVM  3i.iClh<732)I.NCR&'-iENTAL  C=>M  Dje  TO  VANE  MISALIGNMENT 

DCRHAE  ■5I^(TT3  CT5r)TNC^E^ENTAl  "C®!*"  D'JF^TC  5TVIDZL  AS .  "EFFECTS 
DCfMD?  3I3C1P(74u)INCRcM£NTAL  CRM  DJE  TO  DELTA®  VANE  DEFL. 

DCRMD  3“  jTGCmTCT4T)TNryETTTrAL"‘  CRM  DUE  TO  ‘  DELTXT3  VANS  UEFU. 

D  C  9  M  D  R  3I301R(742)INCRSMSNTAL  CRM  DUE  TO  DELTA®  VANE  DEFL. 

D  CRMVH  3IGCTR  C74T5  INCHTMSNTAL  C P*  D~’JE  TO  VTVE  MTS  AL TG NFZNT 

DOME  3I3Q1R(744)INCRE MENTAL  CY  DUE  TC  AEROELAS.  EFFECTS 
iCYDP  "TIGiTTF  C7  45  IT  VCR  EMZ'NTAL  "CT  D"US~TC  DELTA®'  VANE"~trEFTT" 

DC  Y  DU  3IGC1R(7l6)INCSEMENTAL  CY  DUE  'rC  DELTAS  VANE  DEFL. 

DCYDR  31 3CTFT?47>  TNCTcWcMTAL'  CY  DUE  TC  DTUT FT VA N  E"  DEFL." 

OCYMAS  3IG01R(74£> INCREMENTAL  CYM  DUE  TO  AERCELAS.  EFFECTS 
JCYMDP  3IGOTR  (74?  >  I7JT"3  E MENT  AL  CYM  DUE  TO  DELTAF  VANE"  DEFL. 

DCY"<ca  3IS01R  <750  INCREMENT  AL  CYM  DUE  TO  DELTAS  VANE  DEFL. 

OTTTDR  31  j01R(7ST) INCREMENTAL  CTM_DUc~TO  DELTAS  "VANE^TEFUT" 
DCYMVM  3lGC1R<752)lNC5EM£NTAL  CY M  DUE  TO  VANE  MISALIGNMENT 
GAME  3IGG1?Tt757D  EARTH-RELATIVE- FLIGHT  ANGLE  WRT  HORIZONTAL 
G  AM  T  3IG01 R  (752)  ELEVATION  CF  LC  3  TC  TARGET  WRT  LOCAL  NCR. 

3AMX  "  ‘3IGU1  A  (7EV)  ELEVATION"  *rS3ILE  Y-AXTT"  W®T  LCCFL  HOR. 

GRL  3IjC1RC75C)o3CUND  RANGE  LAUNCH  POINT  TO  CURRENT  POS. 

GRT  31  TOT®  (75T  J"S"RC'JNO  '  R'X'NITE  ~FTOW  “C'JjrRFNT  PUSTTTCN  TO  TARGET" 
‘.EDDCT  31301 R  (763)  EARTH-RELATIVE  VELOCITY  IN  LOCAL  NED  FRAME 
5  IG  I  3I"30TS  CTT’T  )  EARTH-PELAT  IVE  FLIGHT  PATH  ANGLE  WFT'NORTH 
EIGT  3I3Q1S <772> AIISUTH  OF  LINE-CF-3I3HT  TC  TARGET 

SIGX  alEOTFtTTT)  ATI‘,TUTH  O"9"  TITS  It  E  X"-XXTT -  -  " 

GRT  3IGC1R<774)SLANT  RANGE  TO  ''HIStT 
U  R  c  «M  31  GO 1R (775) 

VRELM  31 iCI R (73C ) VEL .  OF  MISSILE  RELATIVE  TC  AIR  IN  MISSILE 
XL  oIGGlR (7  34)MI$5ILE  POS.  VECTOP  IN  DESIRED  (T  AR  3  ET) CR  AM  E 

XEDOT  3I3G1R(737)MISSILS  EARTH-RELATIVE  VELOCITY  IN  EARTH  FRAME 
X  7.M  E  3I301P(7PC)P0SITI0N  OF  MIS3ILE  RELATIVE  TC  TARGET  IN  E 
DC  3IGU1°(£1?)FIXED  DSIrTS-X  / Y  t _ A  N  D  Z 

DlV  a  1  G0 1  H  (22 1  )  D  lL T  A  MAGNITUDE  OF  INERTIAL  VELOCITY  . 

OS, MAT  5I301R(?24)3“$2N3ITIVE  DRIFT  PARTIAL  DERIVATIVE  MATRIX 
FJIMS  3IGC1 R (oiS) IMS  FIXED  DRIFT  RATE 

NiF  3IGClR(e3t)IMS^  NEGATIVE  SCALE  FACTORS-X,  Y r  AND  Z 

FSF  "  cIGClR(?39)IMf  POSITIVE"  SCALE  FACTO°S-X*  Y /  AND  Z 

S3DV  2 1 G  0 1  R  (  ?  ■*  2  )  I  m  S  PES£DO»L  DELTA  V  MEASUREMENT 

-.rx  "  il'TOTR  C--.T)  IMS'  AXIS  NCNCRThOGCNALITY-Y  TO  X 

RIX  _  3IG01R  (  !4<i)IMS  AXIS  NCN0RTHC3CNALITY-Z  TO  X 

R «.  Y  at  G0 1  R  (9  47  )  I  *$  AX  I  S  NC  NO  PT  HC  30  N  A  LI  T  Y  -  Z  TO  Y 

TDMG_  31  G0 1  R  (  348  )  M_A£  NI TLDE  CF  IMS  TCRSUcR  RATE 
I XX  Yl301?<%9> 

I  YY  31 3Q1R ('505 

IIZ  '  dI  GvIY  U  31  >  "  "  . 

VCLD  -jIGUI  R (£  52) 

NOZDLX  31 301 R(l:3) NOZZLE  MOMENT  ARM  COMPONENT  ALONG  X  AXIS 
T  C  3  C  N  31301  R(?  SGOTIML  AT  V.MIC*  PCS  JET  WAS  TURNED  ON 
*.,0  71  GC1  R  C  -  *3)-JLZ*  ANGLE  9 »T £3 -* NT £NN A  TO  "ISSILE  FRAME 

'.GO  :lG01t(25e)cULFR  AnK_£  A  C  C 3  . -A  NT £N NA  TO  MISSILE  FRA*: 

TA  CJ  alGOl  R  (  !7L)t:m  E'  RE  iUl;  ED  FOR  CCU  ACQUISITION 


taut*  1-1.  TRW  PROGRAM  OUTPUT  VARIABLES  (Continued) 

TOC  PI  31301  P  (571  >r:y  E  AT  WHICH  C3I  COMMAND  WAS  RECEIVED 
F3IFL3  31 301 1  (  5t‘)  FIRST  STAGE  IGNITION  FLAG 
F3SFL3  si 301  1  (  59>FIRST  STAGE  SEPARATION  FLAG 

101  3IG01K  a?)  0  INPUT  DISCRETE  WORD  1 

102  3 X Q 0 1 1 (  Si)  0  INPUT  DISCRETE  WORD  2 

103  BlSCTlT  69)  O'  INPUT  DISCRETE  WORD  3 

IDEA  3 1  301 1  (  70  )  D -£R  IV  AT IVE  COMPUTATION  COUNTER 

InTTFL  3I31jTIT_71')  INITIALIZATION  FLAG 

ISc3  31301 1  (  725  FLI3HT  SE3LENCE  INDICATOR 

ITTTGE  3 130 1 1  C '777  STA'G'r'PVA" 37=1'  LSuSTCH  AND  STAG E;=CCTH7R WISE 
CD2  31 301  1  (  ?i)  0  OUTPUT  DISCRETE  WORD  2 

rvs f us- Hj orr( t tctpv  Sc?*pArro*TFLAT 

SSIFL5  3I501I (113)SEC0ND  STAGE  IGNITION  FLAG 

S  S  STL'3  1>  1 SOTTTI  T4  J  5  =  C  OW  STAGT- SE9  A’RJTTO'N  T  L  A  S' - 

TiKOPT  31 3 31 1 (11 3)THRUST  0°TI0N  (1-USE  TSLV;  2-USE  TVAC> 

jdflg — snorrczrn  jet  "“flag  tc-no  jd~;  ~i-jd7 - - 

M3T3PT  3IS012C227)J»AS3  TABLE  OPTION  FLAS 

RVJ  “II  GO  TITS'  rZTFrrATTTOr'J  E  f-  GTT7  LAGS  X  2“  j  E  TS7 - - 

A  it  M  C  C  0  31301 1  (435)CrtE  FU  ARM  CODE 

JAC  CET"Tr3UTT(  43  cl  LATTE  RY'  ACTIVATE  TO  57  -  —  -  — -  - 

saFLG  3I3G1I(477)=>ATTcRY  ACTIVATED  FL*3 

cow ro'  JiacncwioysAT  R-TfijiRvo— nrrcF^woRD ~ - - - 

EPwAD  3I301I(439)3AF  REQUIRED  EP  WORD 

ri  R'CGTT  =T  jTTTIT4^CryCTTFFU  ‘TTirT-'C'irDF  .  “ 

FSAAFl  31301  I (441 ) FIRST  STAGE  IGNITION  S  AND  A  APM  FLAG 
F7ACC5 -:»Tu'CTT(4‘4t}  FTR  ST- STAGE'  I'INTTTCN  S  "  AND  A  CTTOf 
IsOFL  3IG31:(445)30C3T  CONFIG.  (1-SINGlEp  2-TWC  STAGE) 

IFAWRT’  sT3CHC444TSAF"RES  J'ITED-  TNFLTGHT "TAMING-  *ORD  . . . 

LOW.RC  31  GG1  1  (445  )  SAF  REQUIRED  LIFT-OFF  WORD 
col  3:T0TT'(w"43>  0  ~  CLTP'JT  DISCRETE 'WORD  1 

P*CNF  51301  I U47JP4C  NO  FAULT  FLAG 
I'VS-mD  323CTlT4'<a)SAF  "ArO'UTPeD'  RV  SEPARATION"  WORD""" 

TII*AD  3IGC1 1  (4a9 J3AF  REQUIRED  TERMINAL  SE3UENCE  INITIATE  WORD 
w.iArtG  7 1  GUI  I  (45G  J  WXfTH7  Ad  A ’MED  FLA'S'  ~  "  “ 

WrtcPLS  313 'J  II  (451  >  W  A  ?  H  s  A  C  EUR  ST  FLAG 

WrtT  jT 3GTIC4  5 2TSXR H Ta'd  "CpTTC n  -  -  - 

CuROV  31 jC1IUo3)CL£AR  IMS  CELT4-V  COUNTERS  f=LAG 

CAY  SC — 31  SOT  IT5"c4Tn‘JM  3ER- j  f  "  0  AY  S  SINCE  LAST  P)S  GAL  16  RATION 

10*.  315012(51 5)  J  IN°UT  DISCRETE  WORD  4 

N3LV  ETGCnTTTS^AJtirATIV'T-pL'LTE  TOjNT-Xp  '  7 /“AT3&- 2 . . . 

OOh  3IGC1 I (532)  0  OLTPUT  DISCRETE  WORD  A 
POLV  3  la  Cl  I  (533)  POSITIVE  »„LS2  COJNT-X,  Yp  AND  2 
TR3TST  3lGC1I(o25)IMS  TORCUER  eIT£  TEST  FlAG 

105  '  TI  SC1TC571}  O  '  INPUT  discrete  WORD  T  ' 

106  S I  301 1 ( c  72  )  0  INPUT  DISCRETE  WOPD  6 

mm  2iscii(67T>d:l  error- yooil  tntti ALirATiofr  flag" 

MNFCNT  31 3011 (c77)MAjCR  F R A ME  COUNTER 

Z’T* - TI  JCTTTT? " 5  JJ  Ou  TP'JT'  D  1 3  C R  zT~~KC  A  D~  l  . 

<rJ AFL  31 30 12(70 6)pIX  DATA  READ  =  f  3  *  C  P- F  LA  G 

m.nTCN  S I  GOT  I  (754)  ANTENNA  ON  -L  AG -STA  r ILI 2 e  ANTENNA  RECEIVED 

DCLUN  3I3u1I(7j5)0:u  ON  FLAG 
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TABLE  1-1.  TRW  PROGRAM  OUTPUT  VARIABLES  (Continued) 


FFiYOCN  irSCTI  (73c)  FTP  ST  S TAG E  HYDRAULIC  SYSTEM  ON  FLAG 
N*VFL3  alGUlI  C7-»2>  =  =  GIN  NAVIGATION  'LAG 

03 - TrSDTXC7So)  73 - 3TJTFuT~t.TTC  ATTs’  WC  RD“3 

PwRFlG  3IG01I(7s7)¥ISSILc  PO*=R-CN  FLAG 

RCSON'  3T5DTI(T6S3  3TAimCNrirCNTR"CL  (TN  FLATT . . .  “ 

SUCN  3I3C1I(7:>9)R4C4R  UNIT  ON  FLA 5 

SHYDTTC  3I3riir77T7)3TC0Nir_5TAC-r^Y0'F"ALLrC  SYSTTK  W  FOG  "  ’ 
•VCSGG1  3ICG1  I(?7c)VCS  GAS  GENERATOR  NJM3ER  1  INITIATE  FLAG 

VGS  GCZ~  7TS  Cl  TT77V  TV  CV  GAS  TETTER  A'T  Jnru»  3&rT~M  IT  1ATE- FLAT? - 

WNDOFT  31  GUI  I  (733)  i»IN  D  OPTION;  G-NO  WINDS/  1-NEV/  2-HAV 
uA«RD'~  3TiH31IC777)3ATTERY  'ACTIVATE  'WCRD'  FROMT  *'ACP  “  ' 

COFLG  3I3C1 I (794) S AF  CUTOFF  FLAG 
£  3  F  L  G — 3IGU T 1X7^57  S  AF  EAPTh  PENETRAT'CTR  FL  AG 

F  5A«RD  3ISC1 1  (796)  FI?  ST  STAGE  IGNITION  5  AND  A  WORD  FROM  PACP 
n  I  ATX.  3TGTJ1  TT7^7TFT7  ST- STA'^E  IT NITITT  “4  PM  ED  Tt  A3 
FSI#PD  31301  1  (79c )  FIRST  STAGE  I3NITICN  WORD  FROM  PACP 
FS3AFL  3TJC1 1  (7^9?  TIP  ST“  STAGE  SEPARATION  ’ARMED  FLAG 
r  SS  «RD  3I301I  COO)FIFST  STAGS  SEPARATION  'WORD  FROM  PACP 
I  FA  r  L  3  TTGOTKojI  73  6  F  INFLIGHT  ARMING  FLAG 
LOFLG  3I3G1 I (3 02) 3 AF  LIFT-OFF  FLAG 

N3  =IGX1 1  (77^1  ^j^R  ’3r TITS  '  EXPECTED  IN  CURRENT  S’A'r  wO’RD 

i.SAa  3IGG1I(C:4)NUP3£R  OF  eATTERY  ACTIVATE  PITS  ACCUMULATED 
!< FS AE  31701  I  (775 )  N'JN  SET  OF  FSI  S  'AND  A' TITS  ACCUMULATED 
\FSIE  _T  3C1 1(3)4) NUF  CE3  OF  FIRST  ST AGE  IGNITION  3ITS  ACCUMULATED 
N  r  S  Go  "  ’317X11  ( ;077  N  JP  E  R  Or  'FTP  ST'  STA  Gs'  ST?.'  PITS  ACCUMULATED 
'.SAFE  EI301KE  03).NUP  =  ER  OF  SAr  cITS  ACCUMULATED 

NS  SIX - JTjUTTTTXTJNUTPER- OT  TTTT_A  CTUTUlATFC  TN”  S  3 ITORD 

NSSGo  ?IjG1I(£1C)NUF3£R  OF  SECOND  STAGE  3E°.  BITS  ACCUMULATED 

SAFRVS  3IGJ1IC31TJGAF  PV  Sc3.  TcA  G  ~  ----- 

SA FmPO  3I3G1I(312)CUFRENT  SaF  WOsD  FROM  PACP 

SSIAFL  31 7CTI(STT) SECOND  STAGE  TGNIT3N  ATME'D  'FLA  3 

aSIwftO  SIGOIKoIOSoCOND  STAGE  IGNITION  WORD  RECEIVED  FROM  ’AC P 

S STTFX -  1T701 1  (  r  11  > 3 'OVTIT A G E  SSPATATICN  ARTD  FLAG 

SaSrfSD  51  G01 i (31 o) SEC  ON  D  STAGE  SEPARATION  WORD  FROM  PACP 

T3IFL3  SIGOlITrrDSAF  TERMINAL  SEQUENCE'  INITIATE  FLAG 

JD23  31 GC 11(354)  0  3REVICUS  VALUE  CF  OUTPUT  DISCRETE  wCRC  2 

acgflg  si gttt  ( Ryu dcl  a  c~cui s i tion  flag  -  ~ 

-05?  31 G01  I  ( -  62 )  0  PREVICU?  VALUE  Or  CJT',UT  DISCRETE  WORD  5 

3XN7L7  JT  3  CTT  ( 739  0  TCATi  '  MED  E"T.A  G  ' 

?  i I o 0 1 C (  54)CISTANC£  FS EARTH'S  CENTER 

3imttm“e: g-ttd c  Fcvsiputation  ttvj  rRo^r  srAir  of ’navigation 

TUcR  31 jC 1 D  (  5E)TIM=  AT  WHICH  DERIVATIVES  ARE  COMMUTED 
X  '37  SC  ID  (TT9)  *TATE  V3GTX'R-Tc7'“’IS  I#  0/  AN  D’  w  VECTORS 

a  3  3IG01C  C143)°3?U0)  (I  NT-E5MEDIATE)  STATE  VECTOR 
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